SPECIFICATION 



TO ALL WHOM IT MAY CONCERN: 

BE IT KNOWN that we, Tetsushi ABE, a subject of Japan 
and residing at Yokohama-shi, Kanagawa, Japan, Hiromasa FUJII a subject of 
Japan and residing at Yokosuka-shi, Kanagawa, Japan, Shigeru TOMISATO a 
subject of Japan and residing at Yokohama-shi, Kanagawa, Japan and Tadashi 
MATSUMOTO, a subject of Japan and residing at Ebina-shi, Kanagawa, 
Japan have invented certain new and useful improvements in 

"TURBO-RECEPTION METHOD AND TURBO-RECEIVER" 
and we do hereby declare that the following is a full, clear and exact 
description of the same; reference being had to the accompanying drawings 
and the numerals of reference marked thereon, which form a part of this 
specification. 
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TURB O-RECEPTION METHOD 
AND TURBO-RECEIVER 



BACKGROUND OF THE INVENTION 
5 The invention relates to turbo-reception method and a 

turbo-receiver as may be used in a mobile communication, for example, and 
which apply an iterative equalization utilizing a turbo-coding technique to 
waveform distortions which result from interferences. 

A task in the mobile station communication business is how to 
1 0 construct a system capable of acquiring a multitude of users on a limited 
: j frequency domain with a high quality A multi-input multi-output (MIMO) 

system is known in the art as means for solving such a task. The architecture 

£ 5 

of this system is shown in Fig. 3 OA where a plurality of transmitters SI to SN 
transmit symbols Ci(i) to c N (i) at the same time and on the same frequency, 
■ 15 and the transmitted signals are received by an MIMO receiver equipped with 
J a plurality of antennas # 1 to #M. The received signals are processed by the 

! : is? 
'. S 

! y receiver, which estimates transmitted symbols ci(i) to c N (i) from the 

respective transmitters SI to SN and delivers them separately as Ci(i) to c N (i) 
to output terminals Out 1 to Out N. 

20 Up to the present time, a study of a specific implementation of an 

MIMO receiver in an MIMO system is not yet satisfactorily warranted. If 
one attempts to construct an MIMO receiver in an MIMO system on the basis 
of MLSE (maximum likelihood estimation) criteria, denoting the number of 
transmitters by N and the number of multi-paths through which a wave 

25 transmitted from each transmitter reaches the MIMO receiver by Q, the 

quantity of calculation required for the MIMO receiver will be on the order of 
2 (Q " 1)N , and will increase even more voluminously with an increase in the 



number of transmitters N and the number of multi-paths Q. When 
information from a single user is transmitted as parallel signals, which are 
then received, a separation of individual parallel signals from each other 
requires a quantity of calculation, which increases exponentially with number 
of multi-paths. Accordingly, the present invention proposes herein an 
improved calculation efficiency turbo-reception method for a plurality of 
channel signals. To start with, an existing turbo-receiver for a single user 
(single transmitter) or a single channel transmitted signal, which illustrates the 
need for the present invention, will be described. 
Turbo-receiver for single user 

An exemplary arrangement for a transmitter and a receiver is 
illustrated in Fig. 31. In a transmitter 10, information series c(i) is encoded 
in encoder 11, and an encoded output is interleaved (or rearranged) by an 
interleaver 12 before it is input to a modulator 13 where it modulates a carrier 
signal, the resulting modulated output being transmitted. The transmitted 
signal is received by a receiver 20 through a transmission path (each channel 
of multipath). In the receiver 20, a soft input soft output (SISO: single-input 
single output) equalizer 21 performs an equalization of delayed waves. At 
the input to the equalizer 21, the received signal is generally converted into a 
baseband, and the received baseband signal is sampled with a frequency 
which is equal to or greater than the frequency of symbol signals of 
information series in the transmitted signal to be converted into a digital 
signal, which is then input to the equalizer 21 . 

For the single user, this corresponds to N=l in Fig. 3 OA, and a 
received output form each reception antenna #m in (m=l, 2, — , M) can be 
represented as follows: 



r m (k) =E q = 0 Q - 1 h m (q) -b(k-q)+v m (k) (1) 
where m represents an antenna index, h a channel value (a transmission path 
impulse response: a transmission path characteristic), b(k-q) a transmitted 
symbol from a user (transmitter 1), and v m (k) an internal thermal noise of the 
5 receiver 20. All outputs from the antennas #1 to #M are denoted by a matrix 
as indicated by an equation (2) to define an equation (3). 



10 where 



r(k)=[r 1 (k)r 2 (k)...r M (k)] T 
=E q = 0 Q - 1 H(q)-b(k-q)+v(k) 

v(k)=[v I (k)v 2 (k)...v M (k)] T 
H(q)=[h 1 (q)...h M (q)] T 



(2) 



(3) 



(4) 
(5) 



ni 



!>!§ 



It is to be noted that [ ] T represents an inverted matrix. In cosideration of the 
15 number of channels Q of the mutipath, the following matrixes and matrix are 
defined: 



20 where 



y(k)=[r T (k+Q-l)r T (k+Q-2)...r T (k) ] T 
=H-b(k)+n(k) 



H = 



'(0) 

0 



H 



'(Q - 1) 
H (0) 



0 



H 



(Q-l)J 



b(k-q)=[b(k+Q-l)b(k+Q-2). . .b(k-Q+l)] T 
n(k) =[v T (k+Q-l)v T (k+Q-2)...v T (k) ] T 



(6) 
(7) 



(8) 

(9) 
(10) 
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r(k) as defined above is input to the SISO equalizer 21 , which is a 
linear equalizer, deriving a log-likelihood ratio Ai(LLR) of a probability that 
each encoded bit {b(i)} is equal to +1 to a probability that it is -1 as an 
equalization output. 



S Pr[b(k) = -l|y(k)] 



^[b(k)] + X p 2 [b(k)] (12) 



sj where A,i[b(k)] represents an extrinsic information fed to a succeeding decoder 

m P 

„ ; 24 and X 2 [b(k)] a priori information applied to the equalizer 2 1 . The 

j 10 log-likelihood ratio Ai[b(k)] is fed to a subtractor 22 where the a priori 

information A, 2 [b(k)] is subtracted therefrom. The result is then fed through a 
deinterleaver 23 to an SISO channel decoder 24, which calculates a 
log-likelihood ratio A 2 as follows: 



! ! i3 



21 S Pr[b(i) = -l|X,[b(i)],i = 0,-,B-l] 

15 B : frame length 

-X 2 [b(i)] + Xf[b0)] (14) 



where X,2[h(i)] represents an extrinsic information which is applied as 

[b(k)] to the equalizer 21 during the iteration, while Aa[b(k)] is applied as a 
priori information A,f [b(i)] to the decoder 24. In a subtractor 25, Xi[b(i)] is 
20 subtracted from A 2 [b(i)] 5 and the result is fed through an interleaver 26 to the 
equalizer 21 and the subtractor 22. In this manner , the equalization and the 
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decoding are iterated to achieve an improvement of an error rate. 

To describe the prestage equalizer 21 in detail, the calculation of a 
linear filter response applied to a received matrix y(k) will be described. 
Using the a priori information X \ [b(k)] for the equalizer 21, a soft decision 
5 symbol estimate 

b'(k)=tan h[X\ [b(k)]/2] (15) 
is calculated. Using the estimate and a channel matrix H, an interference 
component or a replica H*b'(k) of the interference component is reproduced 
and subtracted from the received signal. Thus, 

□ 10 

I y^y^-H-bXk) (16) 

% =H-(b(k)-b*(k))+n(k) (17) 

'■■4 where, 

□ b'(k)=[b , (k+Q-l)...0...b'(k-Q+l)] T (18) 

Because the replica H*b'(k) of the interference component cannot be always a 
correct replica, the interference component cannot be completely eliminated 
by the equation (16). So a linear filter coefficient w(k) which eliminates any 
residue of the interference component is determined according to the MMSE 

20 (minimum mean square error) technique indicated below, 
w(k)=arg min ||w H (k>y'(k)-b(k)|| 2 (19) 
where H represents a conjugate transposition and |j ||a norm. w(k) which 
mimimizes the equation (19) is determined. 

Deriving w(k) in this manner is described in Daryl Reynolds and 

25 Xiandong Wang, "Low Complexity Turbo-Equalization for Diversity 
Channels" (httpr/ee.tamu.edu/Reynolds/). A major achievement of this 
technique lies in a significant reduction in the quantity of calculation. The 
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quantity of calculation according to a conventional MLSE turbo has been 
proportional to the order of 2 Q_1 while a suppression to the order of Q 3 is 
enabled by this technique. It will be seen that w H (k)-y'(k) represents an 
output from the equalizer 21, and is used to calculate A,i[b(k)], which is then 
5 fed through the deinterleaver 23 to the decoder 24 to be used in the decoding 
calculation. 

For purpose of equalization in the equalizer 21, it is necessary to 
estimate the channel value (transmission path impulse response) h appearing 
in the equation (1). This estimation is hereafter referred to as a channel 

1 0 estimation. The channel estimation takes place by using a received signal of 
a known training series such as a unique word which is transmitted as a leader 
of one frame and a stored training series. A poor accuracy of the channel 
estimation prevents an equalization in the equalizer 21 from occurring in a 
proper manner. The accuracy of the channel estimation can be enhanced by 

15 increasing the proportion which the training series occupies in one frame, but 
this degrades the transmission efficiency of the intended data. Accordingly, 
it is desirable that the accuracy of the channel estimation could be improved 
while reducint the proportion of the training series in one frame. 

This is not limited to a receiver for multiple channel transmitted 

20 signals inclusive of MIMO, but the same is true in the channel estimation of a 
receiver such as RAKE receiver or a receiver using an adaptive array antenna 
where the certainty of a decoded result is improved by an iterative decoding 
process. 

The described turbo-receiver has following restrictions: 
25 Tt is an accommodation for a single user (single transmitter) or only 

one series transmitted signal. 

•A channel value (matrix H) is necessary in reproducing an 



interference component, and this must be estimated in actual implementions. 
An estimation error results in a degradation in the effect of an iterative 
equalization. 

It is an object of the invention to provide compensations for these 
5 two restrictions by providing a turbo-reception method and a receiver therefor 
which allow the receiver mentioned above to be extended to a receiver for a 
plurality of transmitted series signals such as for multiple users or parallel 
transmissions from a single user. 

It is another object of the invention to provide a reception method 
10 and a receiver therefor in which a channel value of a received signal is 

estimated from the received signal and a known signal serving as a reference 

;B 

::: signal, the received signal is processed using the estimated channel value, and 

I m 

•1 the processed signal is decoded so that the processing using the estimated 

j; 3 channel value and the decoding are iterated upon the same received signal and 

II 15 which allow the channel estimation to be achieved with good accuracy using a 
relatively short known signal. 
SUMMARY OF THE INVENTION 

According to a first aspect of the present invention, there is 
provided a turbo-reception method of receiving N series (where N is an 
20 integer equal to or greater than 2) transmitted signals. The method 

comprises calculating a channel value h^q) (n= 1, — , N) from M received 
signals r m (m= 1, — , M) and N series known signals, determining a soft 
decision transmitted symbol b' n (k) on the basis of N series a priori 
information ?i 2 [b n (k)] which are obtained by the decoding, and using the 
25 channel value h mrl (q) and the soft decision transmitted symbol b' n (k) to 
calculate an interference component H-B'(k) formed by an intersymbol 
interference produced by an n-th transmitted signal itself and transmitted 
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signals other than the n-th transmitted signal as follows: 



H 



H 



! (0) 



H, 



(Q-l) 



H 



l (0) 



H 



(Q-l). 



H (q) = 



h ll(q) "• h lN(q) 
h Ml (q) ••• h MN(q) 



B'(k)=[b' T (k+Q-l)...b' T (k)...b' T (k-Q+1)] T 
b'(k+q)=[b', (k+q)b' 2 (k+q). . .b* N (k+q)] T 

q=Q-l -Q+l for q^O 

b'(k)=[b , 1 (k). . .0. . .b' N (k)] T q=0 



10 where b'(k) has a zero element at n-th position, Q represents the number of 
multipaths of each transmitted signal wave, q=0, — , Q-l and [ ] T represents a 
transposed matrix, subtracting the intersymbol interference H-B'(k) from the 
received matrix y(k) to obtain a difference matrix y'(k) where 



1 5 y(k)=[r T (k+Q-l)r T (k+Q-2). . .r T (k)] T 
r(k)=[ ri (k)r 2 (k)...rM(k)] T 



determining an adaptive filter coefficient w„(k) applied to the received signal 
for the n-th transmitted signal in order to eliminate residual interference 
20 component in the difference matrix y'(k) using the channel matrix H or 

reference signal, filtering the difference matrix y'(k) according to the adaptive 
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filter coefficient w n (k) to obtain a log-likelihood ratio for the n-th series as the 
received signal for the n-th series transmitted signal from which the 
interference has been eliminated, and decoding by using the the log-likelihood 
ratios for N series. 

5 According to a second aspect of the present invention, in the 

arrangement according to the first aspect fo the invention, when q=0, 

b'(k)=[b'i (k). . .-f(b' n (k)). . .b' N (k)] T 
where an element f(b' n (k)) is located at n-th position, and f( ) represents a 
function which satisfies f(0) = 0 and having as a variable b' n (k) which 
10 satisfies d{f(b\ (k))}/d{b' n (k)}> 0. 

According to a third aspect of the present invention, the 
equalization takes place in a plurality of stages, and the number of series in 
the equalization output is sequentially reduced in the later stages. 

According to a fourth aspect of the present invention, in a 
15 turbo-reception method in which a channel value of a received signal is 

estimated from the received signal and a known signal serving as a reference 
signal, the received signal is processed using the estimated channel value, 
performing a decoding processing upon the processed signal and the 
processing which uses the estimated channel value and the decoding 
20 processing are iterated upon the same received signal, the certainty that the 
decoded hard decision information signal has is determined on the basis of the 
value of an associated soft decision information signal, and a hard decision 
information signal having a certainty which is equal to or greater than a given 
value is also used as a reference signal in the channel estimation of the next 
25 iteration. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram showing a functional arrangement of a system 
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including an embodiment of a turbo-receiver according to a first aspect of the 
present invention; 

Fig. 2 is a schematic view showing a specific example of a 
functional arrangement of a multiple output equalizer 3 1 shown in Fig. 1; 

Fig. 3 is a flow chart of an embodiment of a turbo-reception method 
according to the first aspect of the present invention; 

Fig. 4A is a diagram of an exemplary frame; 

Fig. 4B is an illustration of processings which are performed during 
each iteration, presented in order to illustrate an iterative channel estimation 
method according to a fourth aspect of the present invention. 

Fig. 5 is a diagram of an exemplary functional arrangement for 
deriving a hard decision symbol which is likely to be certain; 

Fig. 6 is a flow chart of an exemplary processing procedure for 
channel estimation according to the present invention; 

Fig. 7A is a diagram of an exemplary functional arrangement of part 
of the equalizer 3 1 according to the second aspect of the invention which 
reflects an error correcting decoding result of a signal being detected; 

Fig. 7B is a diagram showing an exemplary processing procedure 

therefor; 

Fig. 8 is a diagram of an example of a receiver which uses a 
turbo-equalizer iteratively; 

Fig. 9 is a diagram of an exemplary receiver which performs RAKE 
reception-turbo-decoding iteratively; 

Fig. 10 is a diagram of an exemplary receiver which performs an 
adaptive array antenna reception-turbo-decoding iteratively; 

Fig. 11A is a schematic view of a turbo-equalizer; 

Fig. 1 IB is a schematic view of a turbo-decoder; 



Fig. 12 is a schematic view of a receiver which iterates a processing 
of a received signal using an estimated channel and a decoding processing of 
the processed signal; 

Fig. 13 is a flow chart of an exemplary processing procedure of a 
reception method which iterates a processing of a received signal using an 
estimated channel and a decoding processing of the processed signal; 

Fig. 14A is a diagram of an exemplary frame arrangement; 

Fig. 14B is a schematic illustration of an iterative processing of an 
estimation of a channel H and a noise covariance matrix U when a received 
signal contains noises other than white Gaussian noise; 

Fig. 15 is a schematic view of an exemplary functional arrangement 
of part of an equalizer which is used in the estimation of the noise covariance 
matrix U; 

Fig. 16 is a flow chart of an exemplary processing procedure which 
iterates the estimation of a channel value used in the estimation of the noise 
covariance matrix U and the decoding processing; 

Fig. 17 is a diagram showing the principle of a turbo-receiver 
according to a third aspect of the present invention; 

Fig. 18 is a diagram of an exemplary functional arrangement of a 
turbo-receiver according to the third aspect of the present invention; 

Fig. 19 is a diagram showing a specific arrangement of a functional 
arrangement of a multi-user (prestage) equalizer 71 shown in Fig. 18; 

Fig. 20 is a flow chart of an exemplary processing procedure of a 
turbo-reception method according to the third aspect of the present invention; 

Fig. 21 is a diagram showing another exemplary functional 
arrangement of a multi-stage equalization section according to the third aspect 
of the present invention; 
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Fig. 22 is a diagram of an example of a system arrangement in 
which an embodiment according to the first aspect of the present invention (2) 
is applied; 

Fig. 23 graphically shows an error rate characteristic of a 
turbo-receiver according to the first aspect of the present invention (1), with 
assumption that the channel has been perfectly estimated, E b representing a bit 
power and N 0 a noise power; 

Fig. 24 graphically shows an error rate characteristic when the 
iterative channel estimation has been performed while changing a threshold 
value (Th); 

Fig. 25 graphically shows an error rate characteristic of a 
turbo-receiver according to the fourth aspect of the present invention, in 
particular, employing the iterative channel estimation; 

Fig. 26 graphically shows an error rate characteristic of the 
turbo-receiver which employs an estimation of a noise covariance matrix U; 

Fig. 27 graphically shows an error rate characteristic of the 
turbo-receiver shown in Fig. 1; 

Fig. 28A graphically shows an error rate characteristic as plotted 
against E b / N 0 of an embodiment according to the second aspect of the present 
invention which reflects an error correcting decoding result of a signal being 
detected; 

Fig. 28B graphically shows an error rate characteristic plotted 
against a of the embodiment according to the second aspect of the present 
invention; 

Fig. 29 graphically shows a result of a simulation of an error rate 
characteristic of a turbo-receiver according to the third aspect of the present 
invention; 
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Fig. 3 OA is a diagram illustrating the concept of a MIMO system; 

Fig. 3 OB is a diagram of an arrangement in which received signals 
from a pair of antennas are input to a turbo-receiver as four received signal 
series; and 

Fig. 31 is a diagram of a functional arrangement of a conventional 
turbo-transmitter and receiver for a single user. 
EMBODIMENTS OF THE INVENTION 
First aspect of the invention (I) 

Fig. 1 shows an exemplary arrangement of an MIMO system to 
which the present invention is applied. 

In each of N transmitters SI — SN, information series Ci(i) to c N (i) 
are encoded in encoders 11-1,— ,11-N, and the encoded outputs are fed 
through interleavers 12-1,—, 12-N to modulators 13-1,—, 1 3 -N as modulation 
signals, thus modulating a carrier signal in accordance with these modulation 
signals to transmit signals bi(k) to bN(k). In this manner, transmitted signals 
bi(k)-"bN(k) from the transmitters SI, — ,SN form N series transmitted signals. 

A received signal r(k) which is received by a multiple output 
receiver through transmission paths (channels) is input to a multiple output 
equalizer 31. A signal received by the receiver is converted into a baseband 
signal, which is then sampled at one-half the symbol period ,for example, to 
be converted into a digital signal, which is then input to the equalizer 3 1 . 
It is assumed that there are one or more digital signals, the number of which is 
represented by an integer M. For example, received signals from M 
antennas are formed into received signals in the form of M digital signals. 

The equalizer 3 1 delivers N log-likelihood ratios Ai [bi(k)], — 
Ai [t>N(k)]. In subtractors 22-1,—, 22-N, a priori information X\ [bi(k)], — 
X\ [b N (k)] are subtracted from A } [b x (k)] 5 — A] [b^k)], and results are fed 
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through deinterleavers 23-1,— , 23 -N to be input to soft input soft output 
(SISO) decoders (channel decoders) 24-1,— , 24-N to be decoded, whereby the 
decoders 24-1,—, 24-N deliver decoded information series c'i(i)-x' N (i) 
together with A 2 [bi(i)], — ,A 2 [b N (i)] . In subtractors 25-1—, 25-N, Xi [bi(i)], 
— X\ [b N (i)] are subtracted from A 2 [bi(i)], *" ,A 2 [^n(0] > respectively, and 
results are fed through interleavers 26-1—, 26-N to be fed as X 2 [t>i(k)], — X 2 
[b N (k)] to the multiple output equalizer 3 1 and the subtractors 22-1, — , 22-N, 
respectively. 

A received signal r m (k) (m=l,— ,M) from multiple users (a plurality 
of transmitters) is a sum of received signals from the plurality of users when it 
is input to the equalizer 3 1 as indicated below. 



where q=0,— ,,Q-1, and Q represents the number of multipaths for each 
transmitted wave. Defining a matrix y(k) by a similar procedure as used for 
the single user, we have 



r m (k)=Z q=0 Q " 1 £ n =i N Wq)** (k-q)+v m (k) 



(20) 



y(k)=[r T (k+Q-l)r T (k+Q-2). . .r T (k)] T 
=H^B(k)+n(k) 



(21) 
(22) 



where 



r(k)=[ ri (k)...r M (k)] T 



H (0) ••• H(Q-l) 



0 



H = 

o H( 0 ) 



H, 



'(Q-l) 



(23) 



where 



h 



11 (q) h lN(q) 




(24) 



h Ml(q) ••• h MN(q) 



-15- 



R(k)=[b T (k+Q-l). . .b T (k). . .b T (k-Q+l)] T (25) 
b(k+q)=[b! (k+q)b 2 (k+q). . .1* (k+q)] T 
q=Q-l,Q-2,...,-Q+l (26) 

In an interference elimination step, it is assumed that it is desirable 
to obtain a signal from an n-th user (transmitter). In this example, a soft 
decision symbol estimate for signals from all users (transmitters) and a 
channel matrix (transmission path impulse response matrix) H are used to 
produce a synthesis of interferences by user signals other than from n-th user 
and interferences created by the n-th user signal itself, or an interference 
replica H-B'(k) is reproduced. Then , the interference replica is subtracted 
from y(k) to produce a difference matrix y'(k). 



where 



y'(k)=y(k)-H-B«(k) (27) 
=H-(B(k)-B*(k))+n(k) (28) 

B'(k)=[b' T (k+Q-l). . .b ,T (k). . .b ,T (k-Q+l)] T (29) 
b'(k+q)=[b'! (k+q)b' 2 (k+q). . .b' N (k+q)] T 

:q=Q-l,...,-Q+l,q#0 (30) 
b , (k)=[b , 1 (k)...0...b , N (k)] T : q=0 (31) 



b'(k) has a zero element at n-th position. It is to be understood that b' n (k) 
represents a soft decision transmitted symbol estimate which is obtained by 
calculation of b' n (k)= tan h [k 2 [b n (k)]/2] . The matrix B'(k) represents a 
replica matrix of the interference symbol. 

A filter coefficient w n (k) for the n-th user, which is used to 
eliminate the residue of the interference component, namely, a residual 
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interference based on the imperfectness of the interference component replica 
H-B'(k) and an interference component created by the n-th signal itself, is 
determined according to MMSE (minimum mean square error) criteria as one 
which minimizes the following equation (32): 
5 w n (k)=arg min||w n H (k)-y'(k)-b n (k)(| 2 (32) 

Subsequent operation remains the same as for the single user. 
Specifically, w n (k) which is obtained in this manner is used to calculate 
w n H (k)*y'(k), and a result of calculation is fed through deinterleaver 23-n to be 
input as X x [b n (i)] to the decoder 24-n where a decoding calculation is made. 
10 The described method of applying a filter (linear equalization) 

•] processing upon the received signal r m is repeated for all of users 1 to N. As 

i'fj 

l m a consequence, the number of outputs from the equalizer 3 1 will be equal to N, 

and all these outputs are decoded by respective decoders 24-1, 24-N. 
What has been mentioned above is an extension of a turbo-receiver for single 

[ * 15 user to a receiver for multiple users (MIMO). 

!;£ From the foregoing description, it will be seen that an exemplary 

1 y functional arrangement of the multiple output equalizer 3 1 will be as shown in 

Fig. 2. Specifically, M received signals r m (k) are supplied to a received 
matrix generator 311 where a received matrix y(k) is generated and is fed to 
20 equalizers 312-1 to 3 12-N for each user. The channel matrix H which is 
calculated by a channel estimator 28 is also supplied to the equalizers 312-1 
to 3 12-N. A priori information X 2 [b n (k)] from each channel decoder 24-n is 
supplied to a soft decision symbol estimator 313 where a soft decision 
transmitted symbol estimate b * n (k)= tan h [k 2 [b n (k)]/2] is calculated. All of 
25 equalizers 3 12-1 to 3 12-N have identical functional arrangements and process 
in an identical manner, and accordingly, a typical one of them (312-1) will be 
described. 
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In addition, soft decision transmitted symbol estimates b'i(k) to 
b' N (k) are supplied to an interference replica matrix generator 314-1 where an 
interference replica matrix B'i(k) is generated according to the equations (29) 
to (3 1), then the matrix B' { (k) is subject to a filter processing according to the 
channel matrix H in a filter processor 315-1, and a resulting interference 
replica component H-B'^k) is subtracted from the received matrix yi(k) in a 
difference calculator 3 16-lto produce a difference matrix y'i(k). 

At last the channel matrix H or a reference signal which will be 
described later is input to a filter coefficient estimator 317-1 to determine the 
filter coefficient Wi(k) which is used to eliminate the residue of the 
interference component. In the example shown, the channel matrix H and a 
covariance a of a noise component and soft decision transmitted symbols 
b'i(k) to b' N (k) from the soft decision symbol generator 313-1 are input to the 
filter coefficient estimator 3 17-1 to determine the filter coefficient 
Wi(k) which minimizes the equation (32) according to the minimum mean 
square error criteria. A specific example of determining the filter coefficient 
Wi(k) will be described later. The difference matrix y'i(k) is processed with 
the filter coefficient wi(k) in an adaptive filter processor 318-1, and Ai [bi(k)] 
is delivered as an equalized output for the received signal which corresponds 
to the transmitted signal from user 1 . 

A processing procedure for the multiple input multiple output 
turbo-reception method according to the described embodiment of the present 
invention is shown in Fig. 3. At step SI, a channel value h^q) and a 
covariance a 2 of a noise component are calculated from a received signal r(k) 
and each training signal b n (k). At step S2, the channel matrix H is calculated 
from the channel value h^q). At step S3, the soft decision transmitted 
symbol estimate b f n (k)= tan h [A, 2 [b n (k)]/2] is calculated on the basis of the a 
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priori information X 2 [b n (k)] which is obtained during a previous iteration of 
the turbo-reception processing. 

At step S4, the received signal y (k) is generated from the received 
signal r(k). At step S5, the interference replica matrix B' n (k) is generated 
5 according to the equations (29) to (3 1) using the respective soft decision 
transmitted symbol estimates b ' n (k) . At step S6, the interference 
component replica HB f n (k) for the received signal from the n-th transmitter 
is calculated. At step S7, the interference component replica H-B' n (k) is 
subtracted from the received matrix y (k) to provide the difference matrix y' n 
10 (k). At step S8, using the channel matrix H, the soft decision transmitted 
symbols b'i(k) to b^Ck) and the covariance a 2 of the noise component, the 
filter coefficient w n (k) which is used to eliminate the residue of interference 
remaining in the received signal from the n-th transmitter is determined 
according to the minimum mean square error criteria so as to minimize the 
15 equation (32). 

At step 9, a filter processing according to the filter coefficient w n (k) 
is applied to the difference matrix y ' n (k) to obtain the log-likelihood ratio Ai 
[b n (k)]. At step S10., the a priori information X 2 [b n (k)] is subtracted from 
Ai [b n (k)] and the result is deinterleaved and decoded to deliver the 
20 log-likelihood ratio A 2 [b n (k)]. The steps S4 to S10 are performed either 
simultaneously or sequentially for n=l to N. Subsequently, at step SI 1, an 
examination is made to see whether or not the number of decoding operations, 
namely, the number of turbo-reception processings has reached a given 
number. If the given number is not reached, the operation goes to step S12 
25 where the extrinsic information [b n (k)] is subtracted from the 

log-likelihood ratio A 2 [b n (k)], and its result is interleaved to determine the a 
priori information X 2 [b n (k)], thus returning to step S3. If it is found at step 



ry 
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Sl 1 that the number of decoding operations has reached the given number, a 
resulting decoding output is delivered at step SI 3. 

The channel estimator 28 will now be described. Each received 
signal r m (k) can be represented as follows: 

5 r m (k)=2 q=0 Q - 1 L*-i N Uq)'bn (k-q)+v m (k) (33) 

The channel estimator 28 determines the channel value 
(transmission path impulse response) h mn (q) appearing in the equation (33) 
and the mean power (=cr ) of the noise v m (k). Normally, a unique word 
(training signal) which is known to the receiver is inserted at the beginning of 

10 each frame to be transmitted on the transmitting side, as shown in Fig. 4 A, 
and the receiver utilizes the unique word (known signal) as a training series to 
estimate the channel value h mn (q) using RLS (recursive least square) 
technique. Each of the channel decoder 24-1,—, 24-N delivers +1 if a 
log-likelihood ratio A 2 [bi(i)], — A 2 [b N (i)] is positive, and -1 if the latter is 

1 5 negative as a decoded code signal (or transmitted encoded symbol hard 

decision value) b i(i) , b N (i) , which are fed iteratively through interleavers 
27-1,—, 27-N to the channel estimator 28. The received signal r(k) is input 
to the channel estimator 28 together with a unique word which is fed from a 
unique word memory 29 to serve as a reference signal. On the basis of these 

20 input signals, the channel estimator 28 estimate each hnm(q) according to the 
equation (33) and a 2 value according to the least square technique. This 
estimation can be made in a similar manner as the estimation of an impulse 
response when adaptively equalizing the received signal with an adaptive 
filter by estimating an impulse response of a transmission path. 

25 The use of the training series in this manner is conventional, but in 

order to enhance the net transmission rate, it is necessary to reduce the 
proportion that the unique word occupies in one frame, but this increases the 
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error of the channel estimation. If there is such an error, it results in 
degrading the iterative equalization response mentioned above. According 
to the present invention, an iterative estimation of the channel value in the 
following manner is proposed. 

The concept according to the present invention is illustrated in Fig. 
4B. The purpose of the concept is to estimate the channel value iteratively 
in each stage of the iterated equalization processing of the same received 
signal or during the each iteration of the turbo-reception processing. During 
the first run, for information symbol series which follows the unique word, 
the channel value is estimated using only the unique word as a reference 
signal, and the estimated channel value is used to equalize the received signal 
and to estimate the transmitted symbol. However, before the second 
iteration of the equalization processing, the channel estimation is made using 
the unique word as the reference signal, and the symbol estimate (hard 
decision value) which is obtained during the previous decoding processing is 
also used as the reference signal to perform a channel estimation of the entire 
frame. In this instance, rather than using every hard decision value, only that 
hard decision value or values which have been determined as likely to be 
certain is used as the reference signal. The hard decision is rendered by 
delivering +1 when the logic-likelihood ratio A 2 [b n (i)] from the decoder 24-n 
is positive, and delivering -1 if the later is negative. The greater the absolute 
magnitude of the log-likelihood ratio A 2 [b n (i)], it can be said the greater the 
likelihood that the hard decision value is likely to be certain. For example, 
the certainty that is rendered when the log-likelihood 5 is determined to be +1 
is greater than the certainty when the logic-likelihood 0.3 is determined to be 
+1 . Accordingly, an interative channel estimation method which uses a 
threshold value to designate hard decision value b n (i) which is likely to be 
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certain will be described below. 

Initially, using the log-likelihood ration A 2 [b n (i)] from the decoder 
24-n, a soft decision symbol value b' n (i) is determined as follows: 

b' n (i)=tanh[A 2 [b n (i)]/2] 
This operation is made in order to normalize the logic-likelihood value to 1, 
thus preventing its absolute magnitude from exceeding L Next, a threshold 
between 0 and 1 is provided, and when the absolute magnitude of a soft 
decision value b' n (i) is greater than the threshold value, a corresponding hard 

decision value b n (i) is preserved to be used in the iterative channel 
estimation. For example, if the threshold value is chosen to be 0.9, those of 

hard decision values b n (i) which correspond to soft decision values b ' n (i) 
having absolute magnitudes equal to or greater than 0.9 are selected. It is 

believed that the certainty of the selected hard decision value b n (i) is high in 
view of the fact that the threshold value is as high as 0.9, and accordingly, it is 
believed that the accuracy of the iterative channel estimation which is made 
using such hard decision values can be improved. However, it is also 
considered that because the number of selected symbols is reduced, the 
accuracy of the iterative channel estimation may be influenced and becomes 
degraded. Accordingly, it is necessary that an optimum threshold value be 
choosed between 0 and L To add, if the threshold value is to be 1, if follows 

A 

that there is no hard decision value b n (i) which is selected, stating to the 
effect that no iterative channel estimation takes place. Accordingly, a 
threshold value on the order of 0.2 to 0.8 is chosen, as will be further 
described later. 

Accordingly, those of transmitted symbol estimates (hard decision 

values) b i(i) b N (i) for the information symbol series during the first 
transmission which are determined to be likely to be certain according to the 



threshold value are fed from the outputs of the interleavers 27- 1,— , 27-N to a 
previous symbol memory 32 to be stored therein as a previous transmitted 
symbol estimate. During the second iteration of the equalization and 
decoding processing for the received signal r(k) (it being noted that the 
5 received signal r(k) is stored in a memory), the unique word is initially used 
to make a channel estimation, and those of estimated transmitted symbol hard 

/V A. 

decision estimates b i(i) b N (i) which have been determined to be likely 
to be certain are read from the previous symbol memory 32 to be fed to the 
channel estimator 28 to make a channel estimation, namely, a channel 

10 estimation for the entire frame. A resulting estimate tW^q) and a are used 
to perform an equalization and decoding (the estimation of transmitted 
symbol) with respect to the received signal r(k). At this time, those symbol 
values among the estimated transmitted symbols which had been determined 
to be likely to certain according to the threshold are iused to update the 

15 content of the previous symbol memory 32. Subsequently, an estimation 
using the unique word and an estimation using those of the previously 
estimated transmitted signals which are determined to be likely to be certain 
are used to make a channel estimation of the entire frame during a subsequent 
iteration of the equalization and decoding. The estimated channel is used to 

20 perform the equalization and decoding (the estimation of transmitted symbol) 
and to update the previous symbol memory 32. Alternatively, those of the 

transmitted symbol hard decision values bi(i),— , b N (i) from the decoders 
which are determined to be likely to be certain according to the threshold 
value may be directly stored in the previous symbol memory 32 to update it, 
25 and when the symbol values stored in the previous symbol memory 32 are to 
be utilized, they will be fed through the interleavers 27-1,—, 27-8 to be input 
to the channel estimator 28. 
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By the iterations which proceed in this manner, an error of the 
channel estimation can be reduced, the accuracy of the symbol estimation can 
be improved and the problem of a degraded response in the turbo-equalization 
due to an error of an channel estimation can be improved. 
5 When a channel estimation is made with respect to information 

symbol series using the symbol hard decision values which are likely to be 
certain in the manner mentioned above, a function as indicated in Fig. 5 is 
added to each decoder 24-n. The log-likelihood ratio A 2 [b n (i)] is input to 
a soft decision value estimator 241 to calculate b f n (i) = tan h [A 2 [b n (i)]/2] thus 

j; sS, 

10 estimating a transmitted symbol soft decision value b' n (i). The value b' n (i) is 

ij .j. 

j compared against a threshold value Th from a threshold presetter 243, thus 

m 

delivering 1 when b f n (i) is equal to or greater than Th and delivering 0 when it 
is less than Th. On the other hand, the log-likelihood ratio A 2 [b n (i)] is input 
to a hard decision unit 244, which delivers +1 if A 2 [b n (i)] is positive and 
15 delivers -1 if the latter is negative as a symbol hard decision value b n (i). 
The symbol hard decision value b i(i) is fed through a gate 245 which is 
opened when a corresponding symbol soft decision value is equal to or greater 

than the threshold value, and is passed through the interleaver 27-n shown in 
Fig. 1 to be fed to the previous symbol memory 32, thus updating the stored 
20 symbols. 

A procedure of the channel estimation which also utilizes a symbol 
hard decision value or values which are likely to be certain is shown in Fig. 6. 
Initially, at step SI, a channel estimation is made with a received signal r(k) 
and a unique word. At step S2, an examination is made to see if the 
25 decoding processing is for the first time, and if it is, at step S3, the estimated 
channel value h mn (q)is used to perform the equalization and decoding 
processing or the operations shown at steps S3 to S 10 in Fig. 3. 
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At step S4, a transmitted symbol hard decision is rendered with 
respect to the log-likelihood ratio A 2 [b n (i)] to determine a hard decision value 

b n (i). At step S5, b ! n (i)= tan h [A 2 [b n (i)]/2] is calculated from the 
log-likelihood ratio A 2 [b n (i)] , thus estimating a transmitted symbol soft 
5 decision value b ' n (i). At step S6, by examining whether the symbol soft 
decision value b ' n (i) is equal to or greater than the threshold value Th or not, 

those of the symbol hard decision values b n (i) which are likely to be certain 
are determined. At step S7, the symbol hard decision value or values which 
are likely to be certain are used to update the content of the previous symbol 

10 memory 32. At step S8, an examination is made to see whether or not the 
number of decoding operations has reached a given value, and if not, the 
operation returns to step SI, or more exactly, going through the step S12 
shown in Fig. 3 and returning to the step S 1 shown in Fig. 3. 

If it is found at step S2 that the decoding processing is not for the 

1 5 first time, at step S9, a previous stored symbol, namely, a hard decision 

symbol which is likely to be certain is read from the previous symbol memory 
32, and it is used together with information symbol series of the received 
signal r(k) to make a channel estimation, thus transferring to step S3. 

In the above description, even during a second and a subsequent 

20 iteration, the channel estimation is made from the beginning using the unique 
word as a reference signal, but it is also possible that only the hard decision 
symbol which is likely to be certain may be used as the reference signal 
during a second and a subsequent iteration. In this instance, as indicated in 
broken lines in Fig. 6, an examination is made to see if this processing is for 

25 the first time at step ST , and if it is, the unique word is used as a reference 
signal together with the received signal to estimate the channel value. After 
storing the estimated channel value and parameter values which are used in 
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the estimation in a memory at step S3', the operation transfers to the 
equalization and the decoding processing which take place at step S3. 

If it is found at step S 1 ' that the processing is not for the first time, 
before the channel estimation is made, a previously stored channel estimate 
and various processing parameters are set up at S4', whereupon the operation 
transfers to step S9. 

It is to be noted that a solution of the equation (32) is given as 

follows: 

w n (k)=(HG(k)H H +g 2 I) 1 h (34) 
where I represents a unit matrix, o 2 an internal noise power of a receiver 
(covariance of noise component), a 2 I a covariance matrix of noise 
component and G(k) corresponds a squared error of the channel estimation. 

G(k) =E[(R(k) -B'(k) >(B(k) -B'(k) ) H ] 

=diag[D(k+Q-l),. . .,D(k),. . .,D(k-Q+l)] (35) 



where E [ ] represents a mean, diag a diagonal matrix (having elements which 
are not along the diagonal being zero's). 

D(k+q)=diag[l-b' 2 ! (k+q),. . .,l-b |2 n (k+b),. . .,1-b' 2 N (k+q)] (36) 
q=Q-l,Q-2,„.,-Q+l,q#) 
and when q=0, 

D(k)=diag[l-b' 2 1 (k),. . .,1 ,. . .,1-b' 2 N (k)] (37) 
6 1 ' appearing in the matrix D(k) represents an n-th element (it is assumed that 
an n-th user's transmitted signal is assumed to be a desired signal). 
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H 



1 , (Q-l) • N+n 



h = 



H 



'2,(Q-l)-N+n 



(38) 



H 



M-G, (Q-l) ■ N+n 



Thus, h comprises all the elements in the (Q-l)-N+n-th column of H 
appearing in the equation (23). Input to the filter coefficient estimator 317-1 
of the multiple output equalizer 3 1 as shown in Fig. 2 are the channel matrix 
H and the noise power a 2 which are estimated in the channel estimator 28 and 
the soft decision transmitted symbols b'i(k) to b f N (k) from the soft decision 
symbol generator 3 13-1 to calculate the residual interference elimination filter 
coefficient w n (k) according to the equations (34) to (38). 

It will be seen that the equation (34) requires an inverse matrix 
operation, but the required calculation can be reduced by using the matrix 
inversion lemma for the inverse matrix. Specifically, by approximating each 
b' 2 appearing in the equations (36) and (37) by 1, it follows that 



Thus, D(k) have elements having a value of 1 for those elements 
which are located at n-row and n-column while all other elements are equal to 
0. When the error matrix G(k) of the equation (35) which is determined by 
the equations (39) and (40) is substituted into the equation (34), we have 



D(k+q)=diag[0,. . .,0]=0 (q#>) 
D(k)=diag[0,...,l,...,0] 



(39) 
(40) 



w n (kHh'h H +a 2 I)- 1 «h (41) 
where h is as defined by the equation (38). 

With this approximation, w n (k) does not depend on k, and 
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accordingly, an inverse matrix calculation at every discrete time k can be 
dispensed with, thus reducing the quantity of calculation. 

The matrix inversion lemma for the inverse matrix is applied to the 
equation (41). The lemma states that assuming A and B are (M, M) square 
matrices, C a (M, N) matrix and Da(N,N) square matrix, when A is given by 
A'^B" 1 + CD^C 1 *, the inverse matrix of A is given as follows: 

A'^B-BCCD+C 11 BBC)" l C u B (42) 
Applying the lemma to the inverse matrix operation appearing in the equation 
(41), we have 



h(k>h(k) H +c 2 M^+CD'C* 
h(k)-h(k) H =CD- 1 C H ,a 2 I-B 1 , h(k)=C 
I=D^h(k) H =C H 

15 Using these equations to calculate the equation (42), the inverse matrix 
operation appearing in the equation (41) can be calculated. While the 
equation (42) contains an inverse matrix operation (D+ C H BBC) "\ this 
inverse matrix becomes scalar and can be simply calculated. 

Thus, in this instance, it is reduced to the following form: 

20 w n (k)=l/(a 2 +h H h)h (41-1) 

l/( ) on the right side of this equation is scalar or assumes a constant value, 
which may be chosen to be 1 . Thus, we can put w n (k) =h, indicating that w 
(k) can be determined by only h. As indicated in broken lines in Fig. 2, only 
h which is represented by the equation (38) in the channel matrix H from the 

25 channel estimator 28 may be input to the filter coefficient estimator 317-1. 

Without being limited to the use of the matrix inversion lemma for 
the inverse matrix for the approximation by the equations (39) and (40), the 
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approximation also allows the quantity of calculation for the equation (34) to 
be reduced. In particular, when this approximation is used, and the matrix 
inversion lemma for the inverse matrix is used, the quantity of the calculation 
can be further reduced. In addition, if the covariance matrix of the noise 
component is represented by a 2 I, an approximation : w n (k) =h may be used 
as indicated by the equation (41-1), whereupon it is independent from the 
covariance matrix, allowing a further simplification of the calculation. 
Second aspect of the invention (reflecting an error correction) 

In the equalization processing where H-B'(k) is subtracted from the 
received matrix y(k) represented by the equation (27), an error correction 
decoding result is reflected in a transmitted signal soft decision value for a 
signal other than the signal b n (k) being detected, but an error correction 
decoding result which relates to the signal b n (k) being detected is not. For 
this reason, it is preferred to employ the following processing. 

b ! (k) appearing in the equation (29) or the equation (31) is changed 
as follows; 

b'(kHb' 1 (k)b' 2 (k)... 

b' n _! (k>f(b' n (k)) b' n+1 (k). . .b' N (k)] (43) 

where f(b' n (k)) is an arbitrary function having b ' n (k) as an input. 

When such change is made, it becomes possible to reflect an error 
correction decoding result with respect to the signal b n (k) being detected. 
Thus, rather than using b f n (k)=0, by adding a suitable value depending on 
b'n(k), it is possible to emphasize a signal being detected which is buried in a 
noise or an interfering signal, thus allowing b n (k)to be properly detected. 

Because the sign of b ? n (k) relates to a result of a hard decision 



-29- 



rendered upon a symbol which corresponds to b*„(k), and in view of the fact 
that the greater the absolute magnitude of b* n (k), the higher the reliability of 
the hard decision symbol which corresponds to b' n (k), it is necessary that 
f(b' n (k)) satisfies the following requirements: 

The function f has a value of 0 for b' n (k)=0 or when the reliability of 
the hard decision symbol is equal to 0. Or, 



In addition, the greater the value of b'„(k), the greater the value of 
the function £ Or, 

d{f(b' n (k))}/d{b' n (k)}>0 (45) 
Examples of such f(b* n (k)) include the following: 



For example, when the equation (46) is used and a is chosen to be a constant, 
the equation (43) can be implemented in a simple manner. Here, a must 
satisfy the requirement: 0<a<0.6 . If a is greater than 0.6, BER (error rate) 
will be degraded, preventing a correct decoding result from being obtained. 
It is also contemplated to make a to be variable in accordance with the 
reliability of a decoding result. For example, a may be chosen for each 
iteration of the decoding processing. In this instance, the reliability of the 
decoding result will be normally improved with the number of iterations of 
the decoding processings and accordingly, an increased value may be chosen 
for a depending on the number of iterations of the decoding processing. 
Alternatively, the reliability of an entire decoded frame may be determined 
upon each iteration of the decoding processing, and the value of a may be 



f(0)=0 



(44) 



f(b* n (k))=axb' n (k) 
f(b' n (k))=axb' n (k) 2 



(46) 
(47) 
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chosen in accordance with the reliability thus determined. To determine the 
reliability of the decoded frame, a decoding result may be compared against a 
decoding result which is obtained during an immediately preceding iteration, 
and a number of hard decision symbols which changed from the previous 
decoding operation may be counted, for example. Thus, where there is an 
increased number of changed decision symbols, the reliability may be 
determined to low while when the number of changed hard decision symbols 
is small, the reliability may be determined to be high. 

As b' n (k) is changed in this manner, the equation (35) which is used 
when determining the filter coefficient w n (k) for MMSE (minimum mean 
square error) filter is desirably changed as follows: 

G(k)=E[(B(k)-B'(k))-(B(k)-B'(k)) H ] 

=diag[D(k+Q-l),..., D(k),..., D(k-Q+1)] 

Using the equations (29) and (3 1), it follows that assuming 

B'(k)=[b'(k+Q-l)b'(k+Q-2) ... b'(k) ... b'(k-Q+l)] T 
b^kMb'^b^k) ... -f(b' n (k)) ... b' N (k)] T 
b*(k+q)=[b' 1 (k+q)b' 2 (k+q) . . . -f(b' n (k+q)) . . . b' N (k+q)] T 

:q#), q=Q-l, ...,-Q+l 

D(k) has elements located at n-row and n-column, which are represented as 

follows: 

E[(b n (k)+f(b' n (k)))-(b n (k)+f(b' n (k)))*] 
where ( ) represents a complex conjugate. For BPSK modulation, this 
expression turns into the following expression: 
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E[b n (k) 2 +2b n (k)f(b'n (k))+f(b' n (k)) 2 ] 

=E[b n 2 (k)]+2E[b n (k)f(b' n (k))+E[f(b' n (k) 2 ) 

The first term has a mean value of 1 . When b n (k) is approximated by b'(k), 
5 the equation (37) turns into the following form: 

DCk^diagtl-b'^k) l-b'^k^.l-b'^Ck) l+2E[f(b' n (k)b' n (k)) 
+E[f(b' n (k) 2 ) l-o'V^.-.l-b'^k)] (48) 

1 0 For example, when the equation (46) is chosen for f(b' n (k)), D(k) 

turns into the following form: 

D(k)=diag[l-b ,2 1 (k) l-b' 2 2 (k). . . l-b' 2 n .,(k) 

l+(2a+a 2 )b' 2 n (k) l-b'Wk^.l-b'^k)] (49) 

15 

An exemplary functional arrangement which estimates an adaptive 
filter coefficient w n (k) when reflecting an error correction decoding result into 
a signal being detected is shown in Fig. 7 A where a signal being detected is 
chosen as a transmitted signal bj(k) from a first transmitter. A soft decision 

20 transmitted symbol b'^k) is input to a function calculator 331-1 to calculate a 
function f(b'i (k)). Soft decision transmitted signal b'^k) to b' n (k) from N 
decoders and fi^b'^k)) are input to an error matrix generator 332-1 to calculate 
and generate an error matrix G (k) according to the equations (35), (36) and 
(48). The error matrix G(k), an estimated channel matrix H and noise power 

15 a 2 are input to a filter coefficient generator 333-1 where a calculation of the 
equation (34) is made to estimate an adaptive filter coefficient w n (k). In this 
instance, f(b' n (k)) is also input to an interference replica matrix generator 
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314-1, thus generating an interference replica matrix B'(k) represented by the 
equation (29) according to the equations (30) and (43). The filter coefficient 
w n (k) is used in filtering a difference matrix y'(k) in an adaptive filter 3 18-1 , 
thus yielding a log-likelihood ratio Ai [bi(k)]. It will be noted that in the 
5 filter coefficient estimator 317-1 shown in Fig. 2, the function calculator 
331-1 shown in Fig. 7A is omitted, and only the soft decision transmitted 
symbols b^(k) to b' N (k) are input to the error matrix generator 332-1 in order 
to calculate the equation (34). 

In the flow chart of Fig. 3, the interference replica matrix B f (k) is 
□ 10 generated at step S4, and after processings at steps S5 to S7, the filter 

; ; 1 coefficient w n (k) is determined at step S8. If a calculation of the equation 

i! rl 

: (34) is made during the processing at step S8, the soft decision transmitted 

'•j symbols b'^k) to b f N (k) are used to calculate the equations (35) to (37) to 

f 3 generate an error matrix G(k) at step S8-2, and the error matrix G(k), the 

11 15 estimation channel matrix H and the noise power a are used to calculate the 

I5S. 

equation (34) to determine an adaptive filter coefficient w n (k) at step S8-3, as 
shown in Fig. 7B. 

When it is desired to reflect an error correction decoding result in a 
signal being detected as mentioned previously, a soft decision transmitted 

20 symbol b ' n (k)of a signal which is to be detected may be calculated at step S8-1 
before entering the step S4, and this may be used at step S4 where the 
equation (3 1) is replaced by the equation (43), or in other words, the 
equations (29), (30) and (43) may be used to generate an interference replica 
matrix B' (k), and at step S8-2, the equation (37) may be replaced by the 

25 equation (48). In the event f (b ' n (k)) is chosen to be equal to a b ! n (k) or 
a b'n(k) 2 and when a is chosen to be variable, a may be determined in 
accordance with the number of processing operations or the reliability of the 
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entire decoded frame at step S8-1-1, and l+(2a+a 2 )b* n (k) 2 may be calculated 
and used as f (b'„(k)) at step S8-1-2. 

The technique of reflecting an error correcting result into a signal 
being detected is also applicable to a single user turbo-receiver which has 
been described initially in connection with the prior art. In the technique of 
reflecting an error correcting result into a signal being detected, the 
approximation represented by the equations (39) and (40) can be applied, and 
in this instance, only a matrix h shown by the equation (38) which is supplied 
from the channel estimator 28 may be input to the filter coefficient generator 
333-1, as indicated in broken lines in Fig. 7A. 

In the above description, the adaptive filter coefficient w n (k) is 
determined according to the equation (34) or by using the channel matrix H, 
but the use for the channel matrix H can be dispensed with. Specifically, 
during the initial decoding processing (turbo-reception processing), the error 
matrix G appearing in the equation (34) becomes a unit matrix. Accordingly, 
the difference matrix y'(k) and the training signal either alone or in 
combination with a hard decision transmitted symbol b n (k), preferably 
b n (k) having a high reliability in the sense mentioned above are input to the 
filter coefficient generator 333-1 to calculate the adaptive filter coefficient 
w n (k) in a sequential manner by application of RLS (recursive least square) 
technique. Because the error matrix G depends on a discrete time k, during 
a second and a subsequent iteration of the decoding operation, it becomes 
necessary to update the adaptive filter coefficient w n (k) from symbol to 
symbol, and as mentioned previously, it is preferred to determine the adaptive 
filter coefficient w n (k) by using the channel matrix H. 
Fourth aspect of the invention (channel estimation) 

Using not only known information such as a unique word in the 
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iterative channel estimation, but also using a hard decision value of 
information symbol, in particular one which is likely to be certain as a 
reference signal is applicable not only in the described multiple input multiple 
output turbo-reception method, but also generally to a turbo-reception method 
in which a channel (transmission path) of a received signal is estimated from 
the received signal and the known signal, the estimated channel value is used 
to process the received signal and to decode it, and the decoded signal is used 
in iterating the processing according to the estimated channel value and the 
decoding processing upon the same received signal. 

Fig. 8 shows an example in which a hard decision value of 
information symbol is used in a channel estimation turbo-equalizer 4 1 . The 
turbo-equalizer 41 determines a linear equalization filter coefficient in 
accordance with an estimated channel value. The received signal is 
processed by such a linear equalization filter, the processed signal is decoded, 
and the decoded signal is utilized in the iterative processing of the same 
received signal. A received signal r(k) is input to the turbo-equalizer 4 1 and 
is also fed to a channel estimator 42 where a channel value (transmission path 
characteristic) is estimated on the basis of the received signal r(k) and a 
unique word from a memory 29. The received signal r(k) is subject to an 
equalization processing in the turbo-equalizer 41 in accordance with the 
estimated channel value, and is then subject to a decoding processing, 
whereby decoded data c'(i) and a soft decision value b'(i) are delivered. The 
soft decision value b'(i) is input to a symbol designator 43 . If its soft 
decision value b'(i) has an absolute magnitude which is equal to or greater 
than a threshold, the corresponding hard decision value b(i) is stored in a 
previous symbol memory 32 as one which is likely to be certain (as having a 
high reliability) to update it. During a subsequent iterative reception 
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processing (equalization processing) of the same received signal r(k), the 
channel estimation which takes place in the channel estimator 42 takes place 
by using not only the unique word but also the hard decision value b(i) of 
information symbol which is stored in the previous symbol memory 32. 

The turbo-equalizer 41 may comprise the receiver shown in Fig. 1, 
for example, from which the iterative channel estimator 28, the unique word 
memory 29 and the previous symbol memory 32 are removed. It may 
comprise a receiver shown in Fig. 3 1 . Again, a solution of the equation ( 1 9) 
will be in the following form according to the Wiener solution: 

w(k)=E[y'(k)y ,H (k)]-E[b(k)y(k)] 

=[HA(k)H+o 2 I]-h (50) 

where H is as defined by the equation (8), and 

h=[H(Q-l),...,H(0)] T 
where H( ) is as defined by the equation (5) , a 2 =E[ || v || 2 ] (variance of 
noises), and 

A(k)=diag[l-b' 2 (k+Q-1),. . .,1,. . .,l-b' 2 (k-Q+1)] 
In this manner, also in the receiver shown in Fig. 3 1, the channel 
matrix H( ) is estimated, and the channel matrix H ( ) is used to determine the 
equalization filter coefficient w(k), the received signal is filtered according to 
the filter coefficient w(k), and the processed output is subject a decoding 
processing. Accordingly, by using a hard decision information symbol of a 
high reliability in the channel estimation during the iterative recessive 
processing, there can be obtained a more correct channel estimation. 

Fig. 9 shows an example of a turbo-receiver in which the iterative 
channel estimation method is applied in the iterative reception in which 
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RAKE synthesis takes place. A received signal r(k) is fed to a RAKE 
synthesis processor 45 and a channel estimator 42. During an initial 
reception, a channel value is estimated on the basis of the received signal r(k) 
and a unique word in the channel estimator 42, and a compensation for a 
phase rotation which each symbols has experienced on a transmission path as 
well as RAKE synthesis take place in accordance with the estimated channel 
value in the RAKE synthesis processor 45, or a time diversity processing 
takes place to be delivered to a turbo-decoder 46. The turbo-decoder 46 
delivers decoded data c'(i) and a soft decision value b'(i). The soft decision 
value b'(i) is input to a symbol designator 43, and as in described examples, a 

A 

hard decision value b(i) of an information symbol which is likely to be 
certain is stored in a previous symbol memory 32 by updating it. During a 
second and a subsequent iterative reception processing of RAKE reception 
turbo-decoding, not only the unique word but also the hard decision value of 
information symbol which is obtained during the previous iteration are used 
in the channel estimator 42 for purpose of channel estimation, thereby making 
the channel estimation more accurate to improve the quality. 

Fig. 10 shows an example of a turbo-receiver in which the described 
iterative channel estimation method is applied in the iterative reception using 
an adaptive array antenna. A received signal r (k) is received by an adaptive 
array antenna receiver 47, and is then branched into a channel estimator 42 
where a channel estimation takes place in accordance with the received signal 
in combination with a unique word. The estimated channel value is used to 
determine the weights to be applied to each antenna element or corresponding 
reception paths in an array weight determining unit so that the principal beam 
of the antenna directivity response of the adaptive array antenna receiver 47 is 
directed toward the oncoming direction of an intended wave while null is 
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directed toward the oncoming direction of an interfering wave, and such 
weights are applied to applicable locations. A received output from the 
adaptive array antenna receiver 47 is fed to a turbo-decoder 46 to be decoded, 
whereby decoded data c'(i) and a soft decision value b'(i) are delivered. The 
soft decision value b'(i) is input to a symbol designator 43, and a hard 
decision value which is likely to be certain is stored in a previous symbol 
memory 32 by updating it. During a second and a subsequent iterative 
reception processing of the adaptive array antenna receiver 47-turbo-decoder 
46, the hard decision value of information symbol which is obtained during 
the previous iteration is used in the channel estimation in the channel 
estimator 42 together with the unique word. In this manner, the channel 
estimation takes place more correctly, resulting in a more accurate control of 
the antenna directivity response and improving the quality. 

The turbo-equalizer 41 shown in Fig. 8 is schematically shown in 
Fig. 11 A, and as shown, it comprises a series connection of a soft input soft 
output (SISO) equalizer 41a and an SISO decoder 41b, and the operation 
iterates between the equalizer 4 1 a and the decoder 4 1 b. The turbo-decoder 
46 shown in Figs. 9 and 10 is schematically shown in Fig. 1 IB, and as shown, 
it comprises a series connection of a SISO decoder 46a and an SISO decoder 
46b, and an iterative decoding operation takes place between the decoders 46a 
and 46b. The turbo-decoder 46 shown in Figs. 9 and 1 0 may comprise a 
single SISO decoder. 

The examples shown in Figs. 8 to 10 are collectively illustrated in 
Fig. 12. Thus, a received signal is initially processed in an iterative receiver 
(turbo-receiver) 49 in accordance with a channel value which is estimated by 
a channel estimator 42, the processed signal is decoded, and a result of the 
decoding operation is provided in the form of decoded data (symbol) c'(i) and 
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its soft decision value b'(i). The soft decision value b'(i) is compared against 
a threshold value in a symbol designator 43 to determine whether or not a 
corresponding decoded data c'(i) (symbol hard decision value) is likely to be 
certain. If it is determined to be likely to be certain, the hard decision value 
is stored in a previous symbol memory 32 by updating it, and during a second 
and a subsequent iteration of the processing using the estimated channel value 
and the decoding processing, the symbol hard decision value obtained during 
the previous iteration is used in addition to known information such as unique 
word in the channel estimation which takes place in the channel estimator 42 
to provide a more correct channel estimation. 

Fig. 13 shows an exemplary processing procedure of an iterative 
turbo-reception method which also uses a symbol hard decision value. At 
step SI, a channel value is estimated on the basis of a received signal and a 
known signal. At step 2, an examination is made to see if this is the first one 
of the iterative processing, and if it is, the channel value which is estimated at 
step SI is used to process the received signal at step S3, and then the decoding 
processing takes place to determine a symbol hard decision value and a soft 
decision value. At step S4, a symbol hard decision value which corresponds 
to the symbol soft decision value and which is likely to be certain is extracted, 
and at step S5, the extracted symbol hard decision value is used to update a 
previous symbol hard decision value which is stored in the memory 32. At a 
step S6, an examination is made to see if the number of decoding processings 
has reached a given number, and if not, the operation returns to step SI. 
When it is found at step S2 that the iterative processing is not for the first time, 
a previous symbol hard decision value is read out from the memory 32 at step 
S7, and it is used together with an information symbol of a received signal in 
performing a channel estimation, whereupon the operation transfers to step 
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S3. 

Again, as mentioned previously in connection with steps ST to S4' 
with reference to Fig. 6, the known signal may not be used during a second 
and a subsequent iteration. 

In the example shown in Fig. 10, RAKE synthesis processor 45 may 
be inserted between the adaptive array antenna receiver 47 and the 
turbo-decoder 46, as indicated by broken lines. In this instance, a channel 
estimation which is necessary for the compensation of a phase rotation of 
each symbol and the RAKE synthesis in the RAKE synthesis processor 45 
may be served by the channel estimator 42, or may take place separately. 
Noise other than white Gaussian noise 

In the embodiment of the turbo-reception method (according to the 
first aspect of the invention), the embodiment according to the second aspect 
of the invention which takes an error correction into consideration and the 
embodiment of the turbo-reception method characterized in its channel 
estimation method (according to the fourth aspect of the invention), the 
processing took place on an assumption that the noise is white Gaussian noise. 
Specifically, v m (k) appearing on the right side of the equation (29) indicating 
a received signal r m (k) from each antenna is assumed to be white Gaussian 
noise. What is meant by white Gaussian noise is a signal which follows the 
Gaussian distribution and have statistical features expressed as follows: 

E[v m (k)v m (k- q)] =f 2 f q = <> 

0 for q^O 

where E[ ] represents an expected value and a 2 a variance. White Gaussian 
noise may be exemplified by thermal noise which is generated in an antenna 
element. What is influenced by the assumption of white Gaussian noise is a 2 I 
portion appearing in the equation (34) which determines the filter coefficient 
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w n (k) or the equation (50) which determines the filter coefficient w n (k) . For 
example, w n (k) appearing in the equation (34) is calculated through the 
process indicated below. 

w n (k)=(HG(k)H H +E[n(k>n H (k)])"^ 
=(HG(k)H H +a 2 I)" 1 h 

where v m (k) is calculated as E[n (k)-n H (k)] = g 2 I by the assumption of the 
white Gaussian noise having a variance a 2 . The channel matrix H which is 
estimated by the iterative channel estimator 28 (Fig. 1) or 42 (Fig. 12), o 2 and 
the error matrix G(k) which is calculated on the basis of the a priori 
log-likelihood value are substituted into the equation (34) to calculate the 
filter coefficient w n (k). 

When the noise v m (k) is not white Gaussian noise, 
E[n (k)-n H (k)] =a 2 I does not apply. Accordingly, in order to calculate the 
filter coefficient w n (k), it is necessary to estimate an expected value 
(covariance) matrix E[n (k)-n H (k)] for the noise component by a separated 
method. Such method will now be described. A covariance matrix for the 
noise component will be abbreviated as U= E[n (k)-n H (k)] . 
y(k)=H-B(k) + n(k) in the equation (22) is modified into n(k)=y(k) - H-B(k) 
and is substituted into the covariance matrix U, as indicated below. 

U=E[n(k)-n H (k)] 
=E[(y(k)-H-B(k))-(y(k)-H-B(k)) H ] 

If we can assume that a matrix y(k) is available from a received 
signal, an estimate H of a channel matrix H from the channel estimate and 
B(k) is available from a reference signal, it is possible to estimate the matrix 
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U according to the time average method as follows: 

U=2V 0 Tr (y(k)-H-B(k))-(y(k)-H-B(k)) H (51) 
where Tr represents the number of reference symbols. 

During an iterative channel estimation which takes place in the 
iterative channel estimator 28 or 42, the covariance matrix U is estimated 
using the channel matrix H together with the equation (51). A procedure 
therefor is illustrated in Fig. 14. Fig. 14A shows a unique word and 
information symbol series which occurs in one frame of a received signal, and 
Fig. 14B shows the initial processing as well as a subsequent processing. 
During the initial processing, only the unique word is used as a reference 
signal, initially estimating the channel matrix H. Subsequently, U is 
estimated using the unique word and the channel matrix estimate H according 
to the equation (5 1 ). Using the estimate U and H, the filter coefficient w n (k) 
is calculated as follows: 

w n (k)=(HG(k)H H +U)- 1 h (52) 
and the filter coefficient w n (k) is used to apply a first equalization upon the 
received signal, thus estimating transmitted information symbol. 

During a second iteration, the unique word as well as one of 
information signals (*) estimated during the initial equalization which is 
determined to be likely to be certain according to the threshold value are both 
used as reference signals to estimate H again with the same procedure as used 
during the initial processing, thus estimating U again. As this operation is 
repeated, the channel matrix estimate H becomes more accurate with the 
iteration, and the estimate of U becomes more accurate, thus improving the 
accuracy of the filter coefficient w n (k) to improve the response of the 
equalizer. 

In this manner, a turbo-reception when a noise other than white 
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Gaussian noise is contained in a received signal is made possible. 

A functional arrangement in which a linear equalization is 
performed by estimating a covariance matrix U of a noise contained in a 
received signal is shown in Fig. 15 as an example in which a log-likelihood 
ratio Ai [b^k)] is obtained as an equalization output from a multiple output 
equalizer 3 1 shown in Fig. 2 for a received signal corresponding to a 
transmitted signal from a first transmitter. Parts corresponding to those 
shown in Fig. 2 are designated by like reference numerals in Fig. 15 as used 
in Fig. 2. 

A unique word from a unique word memory 29 or a previous 
symbol hard decision from a previous symbol memory 32 which is likely to 
be certain is input to a reference matrix generator 319, which then generates a 
reference matrix B(k) according to the equations (25) and (26). The 
reference matrix B(k), an estimation channel matrix H from a channel 
estimator 28, and a received matrix y (k) from a received matrix generator 
3 1 1 are supplied to a covariance matrix estimator 321 , which then calculates 
the equation (5 1) to obtain an estimated matrix U for a covariance matrix U. 

Soft decision transmitted symbols bi'(k) ,-,b n '(k) from a soft 
decision symbol generator 313-1 are input to an error matrix generator 322-1, 
where an error matrix Gi(k) corresponding to the square error of the channel 
estimation is generated according to the equations (35), (36) and (37). The 
error matrix Gi(k), the estimated covariance matrix U and the estimatied 
channel matrix H are supplied to a filter estimator 323-1, where the equation 
(52) is calculated to estimate a filter coefficient Wi(k). The filter coefficient 
wi(k) and the difference matrix y'(k) from a difference calculator 3 16-1 are 
fed to an adaptive filter 318-1 where a filter processing W!(k) H y'(k) is applied 
to y'(k), and its result is delivered as a log-likelihood ratio Aj [bi(k)]. 
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When reflecting an error correction recording result into a signal 
being detected, a function calculator 33 1-1 as shown in Fig. 7 A is provided, 
as indicated by broken lines in Fig. 15, to calculate f(b' n (k)), and an 
interference replica matrix generator 314-1 uses the equation (43) rather than 
the equation (3 1), and a error matrix generator 322-1 uses the equation (48) 
rather than the equation (37). 

The procedure shown in Fig. 14B is shown as a flow chart in Fig. 
1 6. Specifically, at step S 1 , a received signal r(k) and a known signal (for 
example, a unique word) are used to estimate a channel matrix H, and at step 
S2, an examination is made to see if this processing is for the first time, and if 
it is, the known signal, the estimated channel matrix H and the received signal 
r(k) are used to calculate the equation (5 1) to determine an estimated 
covariance matrix U at step S3. 

At step S4, the estimated channel matrix H, the estimated 
covariance matrix U and an error matrix G(k) which comprises symbol soft 
decision values are used to calculate the equation (52) to estimate a filter 
coefficient w n (k). 

At step S5, the estimated channel matrix H and the filter coefficient 
w n (k) are used to equalize the received signal or to calculate the equation (27) 
to determine w n H (k)-y'(k) to obtain a log-likelihood ratio Aj [b„(k)], 
subsequently performing a decoding process to estimate a hard decision value 
and a soft decision value of a transmitted symbol. 

Purpose of step S6 is to determine a symbol hard decision value 
which corresponds to a symbol soft decision value which exceeds a threshold 
value and which is likely to be certain (or having a high reliability). This 
symbol hard decision value is used to update a symbol hard decision value 
which is stored in a previous symbol memory 32. Subsequently, an 
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examination is made at step S8 to see if the number of the decoding 
processings has reached a given value, and if not, the operation returns to step 
S 1 . However, if a given number is reached, the processing upon this 
received frame is completed. 

If it is found at step S2 that the iterative processing is not for the 
first time, namely, for a second or a subsequent iteration, a symbol hard 
decision value is read from the previous symbol memory 32 at step S9, and is 
used together with information symbol in the received signal to estimate the 
channel matrix H, subsequently transferring to step S3. 

Again, by changing the steps SI and S2 in the similar manner as 
steps SI' to S4' shown in broken lines in Fig. 6, it is possible to avoid the use 
of the known signal during a second and a subsequent iteration. When it is 
desired to reflect an error correction recording result into a signal being 
detected, a function f(b' n (k)) may be calculated at step S 10, as shown in 
broken lines in Fig. 16, and it may be used to obtain the error matrix G(k). 
In either instance, the hard decision transmit symbol may not be used in the 
estimation of the covariance matrix U. The capability to estimate a 
covariance matrix U of a noise contained in a received signal which is other 
than white Gaussian noise is applicable to a variety of useful purposes, as will 
be described below. 

( 1 ) A reception method for a multiple series transmitted signal 
containing an unknown interfering signal is cited. As shown in Fig. 30A, it 
is assumed that an interfering signal i(k) which is unknown to a turbo-receiver 
such as signals from other cells or zones in a mobile communication, for 
example, to be received by the turbo-receiver in addition to N series 
transmitted signal as may be transmitted from N user transmitters. In this 
instance, the equation (20) can be written into the following form: 
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r m (k)=£ q=0 Q_1 S„=i N hmn(q)-bo (k-q)+i(k)+v m (k) (20)' 
In this model, putting i(k)+v m (k)=v f m (k), we have 

r m (k)=2 q=0 Q " ! 2 n =i N Uq)'bn (k-q)+V m (k) (20)" 
Treating v' m (k) as a noise signal other than white Gaussian noise, H and U are 
estimated in a manner mentioned previously, and w n (k) is estimated, and a 
turbo-reception can be made by iterating an equalization of a received signal 
and an estimation of transmitted symbols. 

(2) In a communication system which employs a transmission/ 
reception separation filter, when an oversampling of a received signal is made 
at a higher rate than a symbol period, there occurs a correlation between noise 
components which are contained in received signals which are sampled at 
individual times, and this prevents the noises in the received signals from 
being treated as white Gaussian noise. In other words, the equation (20) 
does not apply Accordingly, an assumption represented as 

E[n(k)-n H (k)]=a 2 I 
does not hold. A processing upon a received signal which is separated by 
the transmission/reception separation filter may utilize equation (51) to 
determine a covariance matrix U, thereby allowing the received signal to be 
properly processed. 

(3) In the described turbo-reception method, the arrangement is 
such that every multipath component from each transmitter (user) on Q paths 
are synthesized. However, in the event there exists a prolonged delay wave 
on channels (for example, assume that paths include one symbol delay, two 
symbol delay and three symbol delay path and there exist separately a thirty 
symbol delay path: in this instance, the thirty symbol delay path component is 
treat as a prolonged delay wave), it is possible to prevent the prolonged delay 
wave from being synthesized, but to treat it as an unknown interference which 
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can be eliminated by an adaptive filter. When the prolonged delay wave 
component is treated as the interfering signal i(k) in the example according to 
the first aspect of the invention (1), it may be eliminated. 

In the processing of a received signal containing a noise other than 
white Gaussian noise, the estimation of the covariance matrix U is applicable 
to a single user turbo-reception method by allowing it to serve in place of 
a 2 I in the equation (50). In a similar manner, it may be used in a RAKE 
synthesis processing reception illustrated in Fig. 9 or a turbo-reception using 
an adaptive array antenna reception shown in Fig. 10, irrespective of a single 
user or multi-user application, or more generally, in a channel estimation in a 
channel estimator 42 and an estimation of the covariance matrix U during the 
iterative decoding operation shown in Fig. 12. For RAKE reception, only 
the channel estimation may be used. 
Third aspect of the invention (multistage equalization) 

In the forgoing description, received signals r x , — ,r M are equalized 
in a multiple output equalizer 3 1 to determine log-likelihood ratios Ai[b(k)], 

A N [b(k)], but in a modification (2) of the first aspect of the present 
invention, there are provided a plurality of equalizer stages in cascade 
connection in a manner such that the number of outputs is reduced toward a 
later stage equalizer. By way of example, Fig. 17 shows that the equalizer is 
divided in two parts where a prestage equalizer (multi-user equalizer) 71 
cancels an interfering component which is located outside an equalization 
range of a poststage, single user equalizer 2 1 At this end, a preprocessing 
including a soft interference cancel and a linear filtering according to 
MMSE(minimum mean square error) criteria, for example, is made, and 
subsequently, the poststage equalizer 2 V performs an equalization of a single 
user having a number of paths equal to Q. 
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Even when an equalization takes place in a cascade manner and a 
linear filter is used in a prestage processing, it is possible to prevent the 
quantity of calculation from increasing prohibitively. 

An arrangement of a multiple output turbo-receiver according to an 
embodiment which is based on this fundamental concept of the first aspect of 
the invention (2) of the turbo-reception method and an exemplary 
arrangement of MIMO system which incorporates the present invention is 
shown in Fig. 18. It is to be understood that parts corresponding to those 
shown in Fig. 1 are designated by like reference characters as used before, 
without repeating their description. (The same holds true in subsequent 
description.) 

Transmitted signals from each transmitter is received through 
transmission paths (channels) by a turbo-receiver 30. The received signal 
r(k) is input to a multi-user equalizer 71, from which signals from N 
transmitters are delivered in the form of signals Ui(k), — , u N (k) , each of 
which is provided in the form it is free from interferences by signals from 
other transmitters, and channel values a^k), — , a N (k) to be input to single 
user equalizers 21-1, -,21-N. These SISO equalizers 21-1, -,21-N deliver 
log-likelihood ratiosAi [bi(k)], Ai [b^k)]. Subsequent processing 
remains similar to Fig. 1. However, the channel values <*i(k), — , a N (k) 
which are used in the single user equalizers 21-1,— ,21-N are channel values 
which are obtained subsequent to the multi-user equalization, and are distinct 
from the channel matrix H. Accordingly, «i(k), — , a N (k) are referred to as 
post-equalization channel information. 

The operation will now be described more specifically. 

Equations (23) to (26) are defined in a similar manner as described 
above in connection with Fig. 1 in consideration of the number of multipaths 
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(channels) Q. 

The purpose of the poststage equalizers 21-1, — ,21-N shown in Fig. 
18 is to equalize intersymbol interference channels by signal symbols of 
respective users themselves [b n (k)> b n (k-l), — , b n (K-Q+l)] (where n=l, N). 
At this end, the prestage equalizer 71 operates to eliminate interferences other 
than [b n (k), b n (k-l), b n (K-Q+l)] (where n=l, N) within y(k). This 
will be discussed below more quantitatively. 

Initially, using the a priori information X 2 P [b n (k)] (where n=l, — ,N) 
of the equalizer 71 which is fed back from the decoders 24-1, — , 24-N, a soft 
decision transmitted symbol estimate b'(k) is determined according to the 
equation (15). 

This soft decision transmitted symbol b ' n (k) and a channel matrix H 
are used to generate a replica H-B'(k) of an interfering signal, which is then 
subtracted from the received matrix y(k). 

y' n (k)=y(k)-H-B'(k) (27)' 
=H-(B(k)-B'(k))+n(k) (28)' 

where 

B'(k)=[b ,T (k+Q-l). . .b' T (k). . .b ,T (k-Q+l)] T (29)' 

b , (k+q)=[b , 1 (k+q)b' 2 (k+q). . .b'„ (k+q). . .b' N (k+q)] T :q=Q-l . .,1 (53) 
b , (k+q)=[b , 1 (k+q)b' 2 (k+q). . .0. . .b' N (k+q)] T :q=0,. . .,-Q+l (54) 

it being noted that b'(k+q) has a zero element at n-th position. 

The operation of subtracting the interference in this manner will be 
referred to hereafter as a soft interference cancel. Assuming that a replica of 
the interfering signal is generated in an ideal manner, it will be seen that y' n (k) 
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resulting from the subtraction can have only the symbol b n (k) of the n-th user 
and an inter symbol interference component caused by the symbol [b n (k-l), — , 
b n (k-Q+l)] of the n-th user itself which results from putting the n-th element 
of b'(k+q) equal to 0 for q=l, — , -Q+l in the equation (54). 
5 In effect, a contribution from the signal of the n-th user (transmitter) 

to the received matrix r(k) is only that resulting that from the symbol [b n (k), 
b n (k- 1 ), — , b n (k-Q+ 1 )] . However, it will be understood from the definition 
of the received matrix y(k) given by the equation (21) that a contribution from 
the signal of the n-th user (transmitter) within the received matrix y(k) which 

10 results as synthesis of the multipaths will contain, when based on the k-th 
symbol b n (k), intersymbol interference components caused by future symbols 
[b n (k+Q-l), b n (k+Q-2), — , b n (k+l)] . Thus, the interference replica includes 
interference components from the future symbols. In this respect, the 
difference matrix y'(k) defined by the equation (27) ' is distinct from the 

15 difference matrix y'(k) defined by the equation (27). 

Accordingly, a next step in the prestage processing in the equalizer 
71 is to eliminate the residue interference which remains after the soft 
interference cancel, namely, a residual interference component which results 
from an imperfect synthesis of the interference replica H-B'(k) and 

20 interference components between future symbols from y ' n (k) using MMSE 
(minimum mean square error) criteria linear filter. In other words, this 
elimination takes place by an arrangement such that a filtering of y' n (k) using 
the filter characteristic w n (k) as indicated by the equation (55) is equal to a 
sum of the symbols [b n (k), b n (k-l), — , b n (K-Q+l)] each multiplied by channel 

25 values a ln (k), a 2n (k), -%a Qn (k). 

w n H (k>y» n (k) « Z q=0 Q_1 a q (k>b n (k-q) 

=a n H (k)-b n (k) (55) 
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Accordingly, what is required is to calculate the equation (55) by determining 
the filter characteristic w n (k) and the post-equalization channel value (channel 
information) a n (k). The derivation of w n (k) and a n (k) will be described. It 
is to be noted that while the filter characteristic w n (k) is distinct from the filter 
5 coefficient w n (k) given by the equations (32) and (34), similar denotations 
will be used for purpose of convenience. 

Desired solutions are defined as solutions of the following 
optimization problem: 

(w n (k),a n (k))=arg min||w n H (k)-y' n (k)-a n H (k>b n (k)|| 2 (56) 
10 provided ai n (k)=l . 

In other words, w n (k) and a n (k) which minimizes the right side of 
the equation (56) are determined. The constraint requirement <Xi n (k)=l is 
added in order to avoid solutions which result in a n (k)=0 and w n (k)=0. While 
solutions can be obtained under the constraint requirement, 
15 ll«n(k)|| 2 =l 

a solution for a ln (k)=l will be described below herein. For brevity, the 
problem will be redefined. Namely, the right side of the equation (56) is 
defined as m n (k) which is minimized in terms of w and a. 

m n (k)=arg min||m n H (k)-Zn (k)|| 2 (57) 
20 provided m n H (k)*e MQ+i— -1 (which is equivalent to ai n (k) =1) and where 

m n (k)=[w n T (k),-a n (k) T ] T (58) 
Zn(kMy n T (k),b(k) n T ] T (59) 
e M Q + i=[0...1...0] T (60) 

25 

it being understood that e M q+i has "1" element at (MQ+l-th position). 
A solution of the optimization problem is given as follows 
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according to Lagrange's method of indeterminate coefficients described in 
literature [2], S. Haykin, Adaptive Filter Theory, Prentice Hall, pp. 220-227; 
m n (k^-Rzz^-eMQ+i/CeMQ+i 11 -Rzz'-eMQ+i) (61) 

where 

R zz =E[z n (k)-z n H (k)] (62) 
E[A] representing an expected value (a mean value) of A. 



= E 



H-A n (k)-H H +o- 2 l H„ 



H 



n 



I 



(63) 



10 A n (k)=diag[D n (k+Q-l),...,D n (k),...,D n (k-Q+l)] (64) 

where I represents a unit matrix, and a 2 noise power (a covariance of white 
Gaussian noise). 



15 



H n - 



hn(Q-l) 0 0 
h n(Q-2) h n(Q _ 1 ) 0 



'n(O) 



h 



n(l) 



0 
0 

Vq-i). 



(65) 



D n (k+q)=diag[ 1 -b' \ (k+q), . . . , 1 -b' 2 n (k+q), . 



l-b ,2 N k+q]] 



:q=Q+l,...,l 



D n (k+q)=diag[l-b' \ (k+q),. . .,1,. . .,1-b' 2 N (k+q)] 

:q=0,...,-Q+l 



(66) 
(67) 



20 



where diag represents a diagonal matrix (all elements other than those located 
along the diagonal of the matrix being zeros). Thus, if the channel matrix H 
and or are known, m n (k) can be determined according to the equation (61). 
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Accordingly, w n (k) and a n (k) can then be determined according to the 
equation (58). 

Using the filter characteristic w n (k) which is calculated in this 
manner, y r n (k) is filtered according to the following equation; 

u n (k)-w n H (k) T ' n (k) (68) 
where H represents a conjugate transposed matrix. 

These filtered n processed results are fed to corresponding 
equalizers 21-n which follow. In this manner, a received signal u n (k) which 
corresponds to the left side of the equation (1) from the n-th user is obtained, 
Om^k) which corresponds to the channel value h mn (q) on the right side of the 
equation (1) is obtained, and the equation (55) which corresponds to the 
equation (1) is determined. Accordingly, a n (k) is applied as an equalizer 
parameter (channel value) to a succeeding equalizer 21-n. This completes 
the prestage processing by the equalizer 71. 

A processing which takes place in the succeeding equalizer 21-n 
and thereafter will now be described. As mentioned previously, because the 
equation (55) corresponds to the equation (1), the operation which takes place 
in the equalizer 21-n for every user may proceed in the similar manner as that 
of the equalizer 21 shown in Fig. 31, and will not be repeated again since it is 
disclosed in literature [1]. Each equalizer 21-n receives u n (k) , <x n (k) and a 
priori information X 2 [b n (k)] from the decoder 24-n as defined above, and 
calculates and delivers outputs a log-likelihood ratio Ai ( LLR) representing 
the ratio of the probability that each encoded bit is +1 and the probability that 
it is -1, as follows: 
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A I [b(k)] = lQg 



Pr[b n (k) = +l|u n (k),k = 0,- > B] 
Pr[b n (k) = -l|u n (k),k = 0,...,B] 



(69) 



-^l[b n (k)] + ^P[b n (k)] 



(70) 



where Xi [b n (k)] represents an extrinsic information fed to a succeeding 
decoder 24-n, and X 2 ? [b n (k)] a priori information applied to the equalizer 3 1 . 
5 Decoder 24-n calculates the log-likelihood ratio A 2 according to the following 
equation: 



10 where X 2 [b n (i)] represents an extrinsic information applied to the equalizer 71 
and the equalizer 21 during the iteration and the X* [b n (k)] a priori 
information applied to the decoder 24-n. The arrangement shown in Fig. 1 8 
performs an iterated equalization and decoding to improve an error rate. 



15 briefly described with reference to Fig. 19. A received signal from each 
antenna is processed in a receiver 70 as a matrix r(k) ^[r^k) — r M (k)] , which 
is used in a received matrix generator 3 1 1 to generate a received matrix y(k) 
according to the equation (21) which takes the individual multi-paths 
(channels) into consideration. 

20 On the other hand, the received signal r(k) from the receiver 70, and 

a known series signal such as a unique word series used for channel 



A 2 [b n (i)] = log 



Pr[b n (i) = +l|Mb n (i)],i = 0,--,B] 
Pr[b n (i) = -l|X 1 [b n (i)],i = 0,-,B] 



(71) 



^2[b n (i)] + ^f[b n (i)] 



(72) 



A functional arrangement of the multi-user equalizer 71 will be 
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estimation which corresponds to each transmitter and which is fed from a 
unique word memory 29 are input to a channel estimator 28 in order to 
estimate a channel matrix H. 

A priori information X* [b n (i)], A/ [b N (i)] are subtracted from 
5 the log-likelihood ratios A 2 [bi(i)], -,A 2 [^(i)] delivered from the respective 
decoders 24-1, — , 24-N to derive extrinsic information X 2 [bi(k)], — , X 2 
[bN(k)] , which is then input to soft decision symbol estimators 313-1, — , 
313-N where soft decision transmitted symbols b'i(k) , — , b' N (k) are 
calculated according to the equation (15) and are then input to an interference 

10 matrix generator 72. In the interference matrix generator 72, a matrix B f (k) 
of symbol estimates which can be interference signals from other transmitters 
are generated for each n according to the equations (29)', (53) and (54). A 
product of these N matrixes B'(k) and the channel matrix H is generated by 
other interfering signal estimators 73-1, 73-N, respectively, thus 

1 5 determining the replica of interfering components H-B(k). 

These N interfering component replicas H-B(k) are subtracted from 
the received matrix y(k) in subtracters 74-1, — , 74-N, respectively, thus 
providing difference matrixes y f i(k), y' N (k). 

The soft decision transmitted symbols b'i(k) , — ,b N(k) are input to 

20 an error matrix generator 75 where error matrices Ai(k), — , A N (k) are 
generated according to the equations (64), (66) and (67). These error 
matrices, the channel matrix H and the noise power a 2 are input to a filter 
characteristic estimator 76 where the filter characteristic w n (k) and the 
post-equalization channel information a n are estimated according to the 

25 equations (58), (60), (61), (63) and (65). These filter characteristics w 1? 

—,w N and difference matrixes y'i(k) , — , y f N (k) are multiplied together in filter 
processors 77-1, — , 77-N, respectively, or the difference matrixes are filtered, 
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thus determining a component u^k), — , u N (k) of the received signal for the 
symbol [b n (k), b n (k-l), — , b n (K-Q+l)] from each user and for each path, from 
which interferences from other user signals are eliminated. These 
components are fed, together with the post-equalization channel information 
5 ai(k), —, a N (k) which are determined in the filter characteristic estimator 76, 
to the single user equalizers 21-1,— ,21 -N shown in Fig. 18. 

A processing procedure for the turbo-reception method according to 
the first aspect of the invention (2) is shown in Fig. 20 where steps 
corresponding to those shown in the procedure of Fig. 3 are designated by like 

10 step numbers as used before. However, the calculation of the interference 
replica matrix B ' n (k) which takes place at step S4 follows the equations (29)', 
(53) and (54). The step S13 uses the soft decision transmitted symbol b' n (k) 
to generate the error matrix A n (k) according to the equations (64), (66) and 
(67). Step S14 uses the channel matrix H, the noise power a 2 and the error 

1 5 matrix A n (k) to determine the residual interference eliminating filter w n (k) 
and the channel information a n according to the equations (58), (60), (61), 
(63) and (65). Step S15 filters the difference matrix y ' n (k) according to the 
residual interference eliminating filter characteristic w n (k) to determine u n (k). 
At step SI 6, a single user equalization is applied to each filtered result u n (k) 

20 to determine the log-likelihood ratio A n [u n (k)], which is then decoded at step 
S 1 0. In other respects, the procedure is similar to that shown in Fig. 3. 

In the forgoing description, the extent of equalization in the 
poststage equalizer 21-n is defined as a zone for intersymbol interference by 
the symbol [b n (k), b n (k-l), b n (K-Q+l)] (where n=l, N), but such extent 

25 of equalization is adjustable. For example, when Q has a very high value, a 
computational load on the poststage equalizer 21-n will be much greater. In 
this instance, the extent of equalization by the poststage equalizer 21-n is 
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chosen to be Q'<Q while the prestage equalizer 71 may be rearranged so as to 
eleminate intersymbol interferences between the signal of the same user 
outside the zone [b„(k), b n (k-l), b n (K-Q'+l)] (where Q r <Q and n=l, N). 
Such modification will be described later. When the equalization is divided 
5 into the prestage and poststage, a previous symbol memory 32 may be 
provided, as indicated in broken lines in Fig. 19, in association with the 

A 

channel estimator 28 so that the hard decision transmitted symbol b n (k) may 
also be used in the estimation of the channel value, thus permitting the 
accuracy of estimation to be improved. 

10 In the example shown in Fig. 17, the prestage multiple output 

equalizer 71 equalizes transmitted signals in N series so that interferences 
from other series may be separated to provide signals u n of N series and the 
post-equalization channel information a n and subsequently, the signal u n of 
each of N series is processed by the post stage single user equalizer 22-n so as 

15 to eliminate intersymbol interference of the same transmitted signal. In this 
manner, the equalization takes place in two stages which are in cascade 
connection. However, a cascade connection of three or more stages may 
also be used. 

By way of example, Fig. 21 shows that a received signal r m of M 
20 series with respect to a transmitted signal of N series is input to a first stage 
equalizer 81 to provide an equalized signal series eri(k) for 1st to U-th 
transmitted series from which interferencs by (U+l)th to N-th transmitted 
series is eliminated and its associated post-equalization channel information 
ea(k), and an equalized signal series er 2 (k) for (U+l)-th to Nth transmitted 
25 series from which interferences by 1 st to U-th transmitted series are 

eliminated and its associated post -equalization channel information eot 2 (k) 
while a second stage includes equalizers 82-1 and 82-2. er^k) and ea^k) 
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are input to the equalizer 82-1 where they are equalized to provide an 
equalized signal series er 3 (k) for 1st -U r th transmitted series among the 1st to 
Uth transmitted series from which interferences by (Ui+l)-th to U-th 
transmitted signal are eliminated and its associated post-equalization channel 
5 information ea 3 , an equalized signal series er 4 (k) for (Ui+l)-th to U 2 -th 
transmitted series among the 1st to U-th transmitted signal from which 
interferences by the 1-st to Ui-th transmitted series and U 2 -th to U-th 
transmitted signal are eliminated and its associated post-equalization channel 
information ea 4 (k), and an equalized signal series er 5 (k) for (U 2 +l)-th to U-th 

10 transmitted series among 1-st to U-th transmitted series from which 
interferences by 1-st to U 2 -th transmitted series are eliminated and its 
associated post-equalization channel information ea 5 (k). 

Similarly, the equalized signal series er 2 (k) and the channel 
information ea 2 (k) are input to the equalizer 82-2 in the second stage to 

15 provide an equalized signal series er 6 (k) and an associated post-equalization 
channel information ea 6 (k), and an equalized signal series er 7 (k) and its 
associated post-equalization channel information ea 7 (k). When N=5, 
equalizers 83-1 to 83-5 in a third stage represent single user equalizers shown 
in Fig. 18. Alternatively, an input equalized signal to the equalizer 83-3 may 

20 comprise two transmitted signals, and the equalizer 83-3 can eliminate mutual 
interferences between the two transmitted signals to provide a set of two 
equalized signals and their associated post-equalization channel information, 
which are in turn supplied to and equalized by next succeeding single user 
equalizers 84-1 and 84-2. As a further alternative, the equalizer 83-4 may 

25 receive the equalized signal er 6 (k) and the channel information ea 6 (k) to 
eliminate mutual interferences for each of the transmitted signals which 
constituted the equalized signal er 6 (k), such as each of three transmitted 
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signals, each of which may be interfered by two other transmitted signals as 
well as an intersymbol interference of each transmitted signal due to 
multipaths. As another alternative, one or more of the equalizers 82-1 and 
82-2 in the second stage may be arranged so that an equalized signal may be 
simultaneously obtained for each of a plurality of transmitted signals. 

What has been described above can be generalized by stating that 
equalizers in a first stage deliver a plurality of equalized signal series and a set 
of post-equalization channel information, and one or more equalizers may be 
provided in each of one or a plurality of stages which are in cascading 
connection for each equalized signal series and its associated set of 
post-equalization channel information so that an equalized output or a 
log-likelihood ratio Ai [b n (k)] may be delivered finally for each of the 1st to 
N-th transmitted series. 

When the equalization takes place in multiple stages which are in 
cascade connection, it is preferred that the number of paths Q for which 
interferences are to be eliminated be reduced toward a later stage so that the 
quantity of calculation can be reduced. In this instance, an interferring 
component from a path which disappears in a later stage be eliminated in an 
immediately preceding equalizer stage. 

An equalization processing which occurs in the arrangement of 
Fig. 21 when a first stage equalizer 81 deals with N transmitted signals, each 
of which has a number of multi-paths equal to Q to provide a group of 
equalized signal series eri(k) comprising U transmitted signals and their 
associated post-equalization channel information eai(k) and a later stage 
equalizer 82-1 performs an equalization for a number of multipaths equal to 
Q' for each transmitted series will be described. 

In a similar manner as described above in connection with the 
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embodiment shown in Figs. 18 and 19, an interference matrix generator 72 
generates an interference matrix B'(k), but the equations used are changed 
from the equations (53) and (54) to equations (53), (54)' and (73) as follows: 

5 b , (k+q)=[b' 1 (k+q)b' 2 (k+q)„^ :q=Q-l,...l (53) 

b'(k+q)=[0...0 b' u+1 (k+q)...b' N (k+q)] T :q=0,...,-Q'+l (54) f 
b'OcfqHb', (k+q)b' 2 (k+q). . .b' n (k+q) 

...b' N (k+q)] T :q=Q',...,-Q+l (73) 

10 The equation (54)' is intended to provide symbols for a first to U-th 

transmitted series themselves, and to provide an equalization except for 
intersymbol interference caused upon each series by itself and relative to each 
other due to multipaths Q' while the equation (73) is intended to eliminate the 
intersymbol interferences upon the 1st to U-th transmitted series themselves 

1 5 and relative to each other due to (Q'+l)-th to Q-th path in as much as the 
number of multipaths is reduced to Q' in the poststage equalization. 

The interferences matrix B'(k) which is obtained in this manner is 
used to generate an interference signal replica H-B'(k) which is then 
subtracted from a received matrix y(k) as follows: 

20 

y' g (k)^y(k)-H*B'(k) (27)" 
=H-(B(k)-B ! (k))+n(k) (28)" 

This operation of subtracting the interference is referred to hereafter 
25 as a soft interference cancel. Assuming that a replica H-B'(k) of an 

interfering signal is generated in an ideal manner, it will be seen that y ? g (k) 
can only have signal components for symbols [b n (k), b n (k-l), 
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b n (k-Q'+l)] (where n=l to U) for 1st to U-th transmitted symbols. 

A residue of interferences which remain after the soft interference 
cancel is eliminated with a linear filter of a MMSE criteria in a similar 
manner as mentioned previously. In this instance, the equation (55) is 
replaced by an equation (55)' indicated below. 

w g H (k)-y' g (k)^E n=1 u E q=0 ^"Wkybn (k-q)=a g H (k)-b g (k) (55)' 



where 



a g (k)=[a uo (k),. . .,ai, Q ..i(k),. . .,a U;0 (k),. . .,a U Q'-i(k)] T (55-1) 
b g (k)=[bj (k),. . .,bj (k-Q'+l),. . .,bu (k),. . .,bu (k-Q'+l)] T (55-2) 



The derivation of w g (k) and a g (k) takes place in a similar manner as 
described previously to determine w g (k) and a g (k) which minimizes the right 
side of the following equation, which stands for the equation (56): 

(w g (k),a g (k))=arg min| |w g H (k)-y' g (k)-a g H (k)-b g (k)|| 2 (56)' 
1 5 provided a i,o(k)= 1 . 

The constraint requirement is added in order to avoid solutions 
which may result in a g (k)=0 and w g (k)=0. While a constraint requirement 
that||a g (k) || =1 may also be used to provide a solution, in the description to 
follow, the problem is rewritten as follows for ai 5 o(k)=l : 
20 m g (k)=arg min||m g H (k)-z g (k)|| 2 (57)' 

provided m g H (k)-e M Q'+i = -1 
where 

m g (kMw g T (k),-a g T (k)] T (58)' 
25 %(k)=[y g T (k),b(k) g T ] T (59)' 

e MQ ' + i=[0... 1 ...0] T (60)' 
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it being understood that eMQ'+i has "1" element at (MQ'+l)-th position. 

A solution of the optimization problem is given as follows 
according to Lagrange's method of indeterminate coefficients disclosed in 
literature [2]; 

m g Rzz^'Cmq'+i /(e M Q-+i H -Rzz'^eMQ'+i ) (61)' 

where 



R zz =E[Z g (k)-Z^(k)] 
= E 



HA n (k)H H +a 2 l Hf 



H D 



I 



(62) ' 

(63) ' 



A n (k)=diag[D n (k+Q-1),. ,D n (k),. . ,D n (k-Q+1)] (64)' 
hl(Q-l) 0 0 hu(Q-l) 0 o 

hi(Q-2) '•• 0 ••■ h u(Q _ 2 ) "•. 0 

hl(Q-3) ; h 1(Q _i) ••• h u( Q_ 3) : h u( Q_i) 



'l(O) 



'U(0) 



hu(Q'-l) 



D n (k+q)=diag[l-b' \ (k+q),. . .,1-b' 2 n (k+q),. . ., 
l-b ,2 N (k+q)] :q=Q+l,...,l 
D n (k+q)=diag[l,. . .,l,l-b' 2 u+ i(k+q),. . 



(65)' 



l-b' 2 N (k+q)] 



:q=0,...,-Q*+l 



D n (k+q)=diag[ 1 -b' 2 ! (k+q), . . . , 1 -b' 2 n (k+q), . . . , 



l-b' 2 N (k+q)] 



:q=Q',...-Q+l 



(66) 
(67)' 
(74) 



Thus when channel parameters are known, m g (k) can be determined 
according to the equation (6 1 )'. In addition, w g (k) and a g (k) (= eai(k)) can 
be determined according to the equation (58)' . Such calculations can be made 
20 in the filter characteristic estimator (76) shown in Fig. 19, for example, and 
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the following equation is calculated in the filter processor 77-1 for purpose of 
filtering: 

er,(k)=w g H <k)y g <k) 
This equalized output eri(k) and the post-equalization channel information 
5 eai(k)=a g (k) are fed to the poststage equalizer 82-1. 

When there are five transmitted series (users), for example, which 
are divided into a group of three transmitted series (users) and a group of two 
transmitted series (users), in the manner mentioned above, the above 
algorithm is carried out with U= 3 and 2, and the two equalized outputs eri(k), 
10 eai(k) and er 2 (k), ea 2 (k) are input to equalizers which are designed to deal 
with the three transmitted series and the two transmitted series, respectively, 
thus obtaining an equalized output for each transmitted series. 

Reflecting an error correction decoding result for a signal being 
detected into a soft decision transmitted symbol in the manner mentioned 
1 5 above is also applicable to a single user turbo-equalizer receiver shown in 
Fig. 8, RAKE synthesis turbo-receiver shown in Fig. 9, a turbo-receiver 
including an adaptive array antenna receiver shown in Fig. 10, a generalized 
turbo-receiver including a channel estimator 42 shown in Fig. 12. 

In Figs. 13, 14 and 15, a symbol hard decision value which has been 
20 determined as being as likely to be certain is also used as a reference signal in 
the estimation of the channel matrix H and the covariance matrix U during a 
second and a subsequent iteration. However, during a second and a 
subsequent iteration, only the unique word may be used as a reference signal 
to utilize the equation (51) to estimate the covariance matrix U while omitting 
25 the estimation of the channel estimation and the estimation of the covariance 
matrix U which utilize the symbol hard decision value. 
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First aspect of the invention (2) (parallel transmission) 

There is a proposal that information series c(i) from a single user be 
transmitted in a plurality of parallel series in order to achieve a high rate 
transmission with a high frequency utilization efficiency. An embodiment of 
5 a turbo-receiver incorporating the present invention which may be used for 
such transmitted signal will now be described. 

Referring to Fig. 22, where parts corresponding to those shown in 
Fig. 1 are designated by like reference characters as used before, on the 
transmitting side, a modulated output signal b(j) from a modulator 13 is input 

10 to a series-parallel converter 14 where each symbol b(j) is sequentially 
distributed into N series. It is assumed that there are series signal 
bi(k) , — , b^k) , the number N of which is equal to or greater than 2. While 
not shown, these signals are transmitted from N antennas after convertion into 
radio frequency signals. 

15 These N series signals are propagated through channels 

(transmission paths) to be received by the turbo-receiver according to the 
present invention. The receiver has one or more receiving antenna, and the 
received signal is input to a multiple output equalizer 3 1 as a baseband digital 
received signal r m (k) (where m=l, 2, M) including one or more (M) signals. 

20 The received signal r m (k) is generated in a manner as shown in Fig. 30B, for 
example. 

The multiple output equalizer 3 1 is constructed in the same manner 
as shown in Fig. 2, and performs according to a processing procedure as 
shown in Fig. 3. Accordingly, an extrinsic information A,i[b(i)] is subtracted 
25 from a log-likelihood ratio A 2 [b(i)] from a decoder 24 shown in Fig. 22 in a 
subtractor 25, and the subtracted output is interleaved by an interleaver 26 to 
provide a priori information A, 2 [b(j)], which is then converted in a 
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series-parallel converter 15 into N series of a priori information X 2 [bx(k)], 
%2 [b N (k)] to be input to the multiple output equalizer 31. 

Accordingly, N series of received signals are subject to a linear 
equalization in the multiple output equalizer 3 1 in the similar manner as 
5 mentioned previously, delivering N log-likelihood ratio series Ai [bi(k)], — 
Ai [b N (k)], which are then input to a parallel-series converter 16 to be 
converted into a single log-likelihood ratio series Ai[b(j)] to be supplied to a 
subtracter 22. With this arrangement, the input signal format to the multiple 
output equalizer 3 1 is similar to that described in connection with Figs. 1 to 3, 

10 and accordingly, N series log-likelihood ratios Ai [bi(k)], — , Ai [b N (k)] can 
be obtained by the equalization which has been mentioned above with 
reference to Figs. 1 to 3, and it will be readily seen that an iterative decoding 
processing is permitted by the use of the series-parallel converter 15 and the 
parallel-series converter 16. In a manner corresponding to a transmitted 

15 signal from a n-th transmitter in the arrangement of Figs. 1 to 3, an n-th or 
(n-th column) transmitted signal among N parallel transmitted signals will be 
equalized. It will be readily seen that the embodiment described above in 
connection with Figs. 4 to 7 is applicable to the reception of the parallel 
transmission of the N series signals. By processing in a plurality of 

20 equalizer stages which are in cascade connection as illustrated in Figs. 18 to 
21, the reception characteristic can be improved as compared to the 
processing through a single equalizer stage as shown in Figs. 1 to 3 1. 

The turbo-reception method and the turbo-receiver according to the 
present invention are also applicable to the reception of convoluted 

25 code/turbo-code + interleaver + multi-value modulation such as QPSK, 

8PSK,16QAM, 64QAM etc.,, TCM (trellis coded modulation) / turbo TCM. 
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Generation of M received signals 

M received signals r^k), — , r M (k) are derived from M antennas #1, 
#M, but may be derived from a single antenna. Alternatively, M (which 
is greater than L) received signals may be obtained from L (which is an 
5 integer equal to or greater than 2) antennas. While not specifically shown in 
Fig. 1, received signals from antennas #1, — , #M are converted into baseband 
received signals ri r M in a baseband converter and sampled to provide 
digital signals ri(k), r M (k) at discrete time k. 

As shown in Fig. 30b, for example, received signals received by 

10 L=2 antennas #1 and #2 may be converted into baseband signals in baseband 
converters 61-1 and 61-2, respectively, the outputs of which may be sampled 
by A/D converters 64-1, 64-2 and 64-3, 64-4 in accordance with a sampling 
signal from a sampling signal generator 62 and a sampling signal which is 
obtained by shifting the phase of the first mentioned sampling signal by one 

1 5 half period thereof T/2 in a phase shifter 63 to provide digital signals rj(k), 
r 2 (k), r 3 (k), r 4 (k), respectively. Then the digital signal may be input to a 
turbo-receiver 30 as shown in Fig. 1 or Fig. 18 or Fig. 22 to provide N 
decoded outputs. It should be understood that the frequency of the sampling 
signal from the sampling signal generator 62 is chosen so that each sampling 

20 period of the received signals r t (k), • • • , r 4 (k) input to the turbo-receiver 3 0 
coincides with the sampling period which is used when a single received 
signal r M (k) is received per antenna. 
EFFECTS OF THE INVENTION 

As discussed above, according to the first aspect of the invention 

25 (1) , there is realized a multiple output (MIMO) reception method. To 

illustrate a quantitative effect , an error rate response is graphically shown in 
Figs. 23 and 24. In each Figure, E b /N 0 on the abscissa represents a bit power 
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to noise ratio. A simulation presumes following conditions: 



Number of users (transmitters) 2 
The number of multipaths of each user Q 5 
5 Number of reception antennas 2 

Number of information symbols per frame 450 bits 
Number of unique words per frame 25 bits 

Channel estimation method RLS (obliteration factor 

0.99) 

10 Error correction code rate 1/2, constrained length 3 

convoluted code 

Doppler frequency 1000Hz (Rayleigh fading) 

Modulation BPSK 

Transmission rate 20Mbps 
15 Decoder 24 Max-Log-Map decoder 

Number of iterations 4 

No fading within a frame 

No approximation by the matrix inversion lemma for the inverse matrix is 

20 used in the calculation of the filter coefficient w. 

Fig. 23 shows an error rate characteristic when the channel 
estimation is achieve perfectly (without no estimation error) or when the 
channels are known for a number of users N=2, a number of reception 
antennas M=2, a number of Raleigh paths Q=5. It will be seen that initial 

25 run is not an iteration and the second iteration is the first iteration. It will be 
seen that the error rate characteristic is significantly improved by the iteration. 
It will be seen from this that the turbo-reception method according to the 
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invention for MIMO operates in a proper manner. 

Fig. 24 illustrates the effect of an iterative channel estimation 
(according to the fourth aspect of the invention). The abscissa indicates a 
threshold value Th. E b /N 0 is fixed to 4 dB (where E b is for one user), and 
5 Th=l .0 may be considered as a conventional method in which no symbol hard 
decision value is selected, or where a channel estimation using a symbol hard 
decision value is not made. In this instance, it will be seen from the drawing 
that BER characteristic has a little effect of iteration because the channel 
estimation is inaccurate. A threshold Th=0 indicates that all of hard decision 

10 values are directly used, and when hard decision values of information 

symbols are used in this manner, it will be apparent from the drawing that a 
mean bit error rate is improved, allowing the channel estimation to be 
performed exactly to a corresponding degree. For a threshold value on the 
order of 0.2 to 0.6, the mean bit error rate is reduced than when Th=0, 

15 indicating that it is preferable that only a hard decision value which is likely 
to be certain be used. It will be seen that Th which is around 0.25 is most 
preferable. 

Fig. 25 shows an error rate characteristic of an MIMO reception 
method in which a transmitted symbol hard decision value which is 

20 determined to be likely to be certain according to a threshold value is used in 
the channel estimation or which employs an iterative channel estimation in a 
form of a curve 66. In this instance, the threshold is set up at 0.25 and the 
result shown represents a characteristic after four iterations with parameters 
that N=2, M=2, Q=5Rayleigh, f d T s =l/20,000, and 900 symbols/frame. For 

25 purpose of comparison, an error rate characteristic with a perfect channel 

estimation is shown as a curve 67 while an error rate characteristic in which a 
hard decision value of information symbol is not used in a channel estimation 
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or only a single channel estimation is made without iteration is shown by a 
curve 68. It will be seen from this graphic illustration that when the iterative 
channel estimation is used, the error rate characteristic approaches to that 
obtained by a perfect channel estimation. 



determining whether a hard decision value is or is not likely to be certain on 
the basis of a decoded soft decision value, and by using symbol information 
having a hard decision value which is likely to be certain in the channel 
estimation during the next iteration, the channel estimation can be performed 
10 more correctly, allowing a decoding quality to be improved. 

In order to confirm the effect of an embodiment in which a 
covariance matrix U (for noise other than Gaussian noise) is estimated, a 
simulation is made with parameters indicated below. 

1 5 Number of users (transmitters) 3 (one of which is chosen as 



5 



With the channel estimation method mentioned above, by 



an unknown interference: i(k)) 



20 



The number of multipaths of each user Q 5 
Number of reception antennas 3 
Number of information symbols per frame 450 bits 
Error correction code rate 1/2, constrained length 3 



convoluted code 



Doppler frequency 



1000Hz (Rayleigh fading) 



Modulation 



BPSK 



25 



Transmission rate 



Decoder 24 



20Mbps 

Log-MAP decoder 



Number of iterations 



4 



Three users (transmitters) 



are chosen to be of an equal power. 
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Fig. 26 shows a result of simulation of BER (bit error rate) 
characteristic of a turbo-receiver which estimates H and U as shown in Figs. 
14, 15 and 16, and Fig. 27 shows BER characteristic which directly uses the 
turbo-receiver shown in Fig. 1 (which employs the method shown in Fig. 13). 
5 In Fig. 26, it is assumed that the noise comprises only white Gaussian noise, 
and it is seen that two or more iterations of the channel estimation and the 
decoding processing has little effect. By contrast, it will be seen from Fig. 
27 that as the number of iterations is increased, an improvement in the BER 
characteristic is achieved, and in addition, BER is considerably less than that 

10 shown in Fig. 26 for the same E b /N D . 

In order to confirm the effect of the embodiment (according to the 
second aspect of the invention) in which an error correction decoding result is 
reflected into a symbol soft decision value b' n (k) of a received signal from an 
intended user (transmitter), a simulation is made with parameters as indicated 

15 below 

Number of users (transmitters) N 4 
The number of multipaths of each user Q 5 
Number of reception antennas 2 

20 Number of information symbols per frame 900 bits 

Error correction code convoluted code (coded rate: 1 12, 

constrained length 3) 
Modulation BPSK 
Decoder 24 Log-Map decoder 

25 Erroneous coded rate 1/2 

Number of iterations 5 
A choice is made that f(b'„(k))= axb' n (k) 
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Fig. 28 shows BER characteristic of a multiple output 
turbo-receiver shown in Fig. 1 and a multiple input multiple output 
turbo-receiver in which an error correction decoding result is reflected into 
b' n (k), with plotted points shown in black for the former and in white for the 
5 latter. It is to be noted that in the plot, a circle represents an initial run, a 
down-directed triangle a second iteration, a diamond a third iteration, a 
left-directed triangle a fourth iteration and a right-directed triangle a fifth 
iteration. Fig. 28A shows a result of a simulation for the BER characteristic 
plotted against Eb/N 0 when a is fixed to 0.2, and Fig. 28B shows a result of 

1 0 simulation for the BER characteristic plotted against a when Eb/N c is 6dB. It 
will be seen that a=0 represents that b' n (k)=0. It is noted from Fig. 28A that 
in the multiple input multiple output receiver in which an error correction 
decoding result is reflected into b' n (k), an improvement in BER after a third 
and a subsequent iteration is greater in comparison to BER obtained during a 

15 preceding iteration, as compared with the multiple input multiple output 
turbo-receiver shown in Fig. 1, and that after a third and a subsequent 
iteration, when Eb/N 0 required to achieve each BER in a range BER>10" 4 is 
compared, the multiple input multiple output turbo-receiver in which an error 
correction decoding result is reflected into b' n (k) has a gain of about 0.5 dB or 

20 greater in comparison to the multiple input multiple output turbo-receiver 

shown in Fig. 1 . It is also noted that during the fifth iteration at E b /N Q =6dB, 
BER=10" 5 is attained, which is a reduction in BER by a factor of 1/10 or more 
in comparison to the receiver shown in Fig. 1 . It will be seen from Fig. 28B 
that an improvement is achieved in a range of a indicated by an equality 

25 0<a<0.6 and that when a exceeds 0.6, the BER characteristic becomes 

degraded, preventing a correct decoding result from being obtained. From 
this result, it is seen that an optimum value of a is 0.2. However, it should 
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be understood that the value of a is not limited to the optimum value, but a 
suitable range of a having an improving effect can be changed depending on 
the number of users to be received, a propagation environment having 
interferences, a number of receiving antennas or the like. In addition, a 

5 different value may be chosen as an optimum value for a. 

When BPSK modulation is used for a number of users (transmitters) 
which is equal to N, a number of multipaths from each transmitter which is 
equal to Q and a number of receiver antennas which is equal to M, the 
quantity of calculation which is required in an equalizer when a conventional 

1 0 single user turbo-receiver is directly extended to a multiple output (MIMO) is 
on the order of 2 N(Q_1) , as mentioned previously, but with the turbo-reception 
method according to the third aspect of the invention, the quantity of 
calculation can be reduced to the oreder of N(MQ) 3 . By way of example, 
assuming that N=8, Q=20 and M=8, 2 N(Q " 1) ~5-10 45 while N(MQ) 3 ~37-10 7 , thus 

1 5 demonstrating that the quantity of calculation can drastically be reduced 

according to the turbo-reception method according to the second aspect of the 
present invention. 

A simulation has been conducted under conditions given below in 
order to confirm that a good bit error rate characteristic can be obtained 

20 according to the turbo-reception method according to the third aspect of the 
invention. It is assumed that the channel matrix H is known. 



Number of users N 4 

The number of multipaths of each user Q 5 

25 Number of reception antennas 2 

Number of information symbols per frame 900 bits 
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Error correction code 

Doppler Frequency 
Modulation 
Decoder 24 
Transmission rate 
Decoder 

Number of iterations 



rate: 1/2, constrained length 3 
convoluted code 
1000Hz (Rayleigh fading) 
BPSK 

Log-Map decoder 
20Mbps 

Log-MAP decoder 
6 



10 It is assumed that the channel estimation takes place in an ideal 

manner. 

Fig. 29 shows a result of simulation for BER (bit error rate 
characteristic). The abscissa represents mean Eb (bit power)/ N 0 (noise 
power), denotation f d in the graph represents a doppler frequency and D s the 

1 5 period of a transmitted symbol. MRC shown on this graph is BER 
characteristic upon Viterbi decoding of a signal after a maximal ratio 
combining on an order 10 (two antennas><5 paths) diversity channels, and this 
corresponds to the BER characteristic which results when the equalizer has 
perfectly canceled interferences. Thus, the quality of the receiver can be 

20 evaluated by seeing how the BER after the iterations is located close to the 
MRC curve. It is seen from Fig. 29 that with the turbo-reception method 
according to the second aspect of the present invention, the higher E b /N G , the 
less BER, and the greater the number of iterations, the more the BER 
characteristic approaches MRC BRE characteristic, in particular BER 

25 characteristic is located very close to MRC with a number of iterations equal 
to 6. Thus, it is confirmed that the multiple output turbo-receiver which 
employs the turbo-reception method according to the third aspect of the 
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present invention operates properly under severe conditions of four 
each having five paths and using two reception antennas. 



